I  • 


AD-A112  630  SYRACUSE  UNIV  N.Y  DEPT  OF  ELECTRICAL  ENGINEERING  P/G  20/3 

ELECTROMAGNETIC  C0UPLIN6  THROUGH  APERTURES. CU) 

JAN  62  J  R  MAUTZp  R  F  HARRINGTON  F30602-79-C-0011 

RADC-TR-81-380  -NL 


UNCLASSIFIED 


IE  COPY 


kaoc-tr-8  i-aao 

Phase  import 
Jonuory  1982 


ELECTROMAGNETIC  COUPLING 
THROUGH  APERTURES 

Syracuse  University 


Joseph  R.  Moufx 
Roper  F.  Horrlnp*©" 
David  T.  Auckland 


ROME  AIR  DEVELOPMENT  CENTER 
Air  Force  System*  Commend 
Griffis*  Air  Force  Bose,  New  York 


13441 


This  report  has  been  reviewed  by  the  RADC  Public  Affairs  Office  (PA)  and 
is  releasable  to  the  National  Technical  Information  Service  (NTIS).  At  NTIS 
it  will  be  releasable  to  the  general  public,  including  foreign  nations. 

RADC-TR-8 1-380  has  been  reviewed  and  is  approved  for  publication. 


APPROVED: 

ROY  F.  STRATTON 
Project  Engineer 


APPROVED: 

EDMUND  J.  WESTCOTT 
Technical  Director 

Reliability  &  Compatibility  Division 


FOR  THE 


Acting  Chief,  Plans  Office 


If  your  address  has  changed  or  if  you  wish  to  be  removed  from  the  RADC 
mailing  list,  or  if  the  addressee  is  no  longer  employed  by  your  organisation, 
please  notify  RADC.(RBCT)  Griffins  APB  NY  13441.  This  will  assist  .us  in 
maintaining  a  current  mailing  list. 

Do  not  return  copies  of  this  report  unless  contractual  obligations  or  notices 
on  a  specific  document  requires  that  it  be  returned. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  This  pace  (Wh* i  0.1.  Enfrod) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


4.  TITLE  (mnd  SubtUlm) 

ELECTROMAGNETIC  COUPLING  THROUGH  APERTURES 

5.  TYRE  OF  REPORT  *  PERlOO  COVERED 

Phase  Report 

6  PERFORMING  O^G.  REPORT  NUMBER 

N/A 

7.  AuTMORfA) 

Joseph  R.  Mautz 

Roger  F.  Harrington 

David  T.  Auckland 

4.  CONTRACT  OR  GRANT  NUMBERC) 

F30602-79-C-0011 

9.  PERFORMING  ORGANIZATION  NAME  AND  AOORESS 

Syracuse  University 

Department  of  Electrical  And  Computer  Engrg. 
Syracuse  NY  13210 

10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  4  WORK  UNIT  NUMBERS 

62702F 

23380317 

1  ».  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

Rome  Air  Development  Center  (RBCT) 

12.  REPORT  DATE 

January  1982 

Griffiss  AFB  NY  13441 

13.  NUMBER  OF  PAGES 

128 

14.  MONITORING  AGENCY  NAME  4  AOORESSfU  diffarant  from  Controlling  Olll co) 

Same 

15.  SECURITY  CLASS,  (ot  (him  rmport) 

UNCLASSIFIED 

IS..  DECLASSIFICATION'  DOWNGRADING 
N/ASCHE0ULE 

16.  DISTRIBUTION  STATEMENT  ( ol  thlm  Rmport) 

Approved  for  public  release;  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (of  thm  mbmtrmct  mntmrmd  in  Block  20,  It  dlHmrmnt  from  Rmport) 

Same 

16  supplementary  notes 

RADC  Project  Engineer:  Roy  F.  Stratton  (RBCT) 

19.  KEY  WOROS  ( Contlnum  on  rm*mrmm  mtdo  it  n«««n«y  an d  Identify  by  block  r-tnbtr) 

Aperture  Admittance  Ge--pr  ' 

Aperture  and  Conducting  Body  f; 

Aperture  in  a  Screen  Hu.>. 

j^avity-backed  Aperture  Wavegu* 

ized  Network  Parameters 
jf  Moments 

Jots 

ue-fed  Apertures 

'Mi  ABSTRACT  (Canllmi*  on  tomormo  mldo  It  nocmmmofy  md  tdmntUs  by  block  numbor ) 

^This  report  is  a  tutorial  presentation  on  the  use  of  the  generalized  net¬ 
work  approach  to  aperture  problems.  The  general  theory  is  given  and 
applied  to  the  cases  of  electromagnetic  transmission  through  an  aperture 
in  a  conducting  plane,  a  waveguide-fed  aperture,  a  cavity-backed  aperture, 
a  narrow  slot,  and  an  aperture  backed  by  a  conducting  body.  The  last 
two  problems  demonstrate  the  phenomena  of  aperture  resonance  for  slots 
in  thick  conductors  and  apertures  backed  by  conducting  bodies.  At 

DO  i  J»H  71  1^73  EDITION  OF  I  NOV  •»  IS  OBSOLETE 

/ 


UNCLASSIFIED 

SECURITY  CL  AJSlFlCATlON  OP  THIS  PAGE  (Ph-n  Dmtm  Entotod) 


CONTENTS 


Page 

CHAPTER  1  BASIC  CONCEPTS -  1 

1-1  Introduction -  1 

1-2  Fields  and  Sources -  2 

1-3  Potentials -  4 

1-4  Potential  Integrals -  6 

1-5  Conditions  at  Surfaces -  7 

1-6  The  Equivalence  Principle - 10 

1-7  Integral  Equations -  12 

1-8  Method  of  Moments -  15 

1- 9  References  for  Chapter  1 -  18 

CHAPTER  2  A  GENERALIZED  NETWORK  FORMULATION  FOR  APERTURE 

PROBLEMS -  19 

2- 1  Introduction -  19 

2-2  General  Formulation -  20 

2-3  Linear  Measurement -  25 

2-4  Transmitted  Power -  28 

2-5  Discussion -  29 

2- 6  References  for  Chapter  2 -  30 

CHAPTER  3  APERTURES  IN  PLANE  CONDUCTORS -  31 

3- 1  General  Theory -  31 

3-2  Transmission  Parameters -  35 

3-3  Admittance  Matrix -  36 

3-4  The  Rectangular  Aperture -  38 

3-5  Plane  Wave  Excitation  and  Measurement  Vectors -  40 

3-6  Representative  Computations -  41 

3-7  Discussion -  45 

3-8  References  for  Chapter  3 -  48 

iii 


Page 


CHAPTER  4  WAV ECU I DE-FED  APERTURES -  48 

4-1  General  Theory -  48 

4-2  Application  to  a  Rectangular  Waveguide -  55 

4-3  Sample  Computations -  57 

4- 4  References  for  Chapter  4 -  60 

CHAPTER  5  CAVITY-BACKED  APERTURES -  61 

5- 1  Theory  for  a  Cylindrical  Cavity -  61 

5-2  Aperture  Backed  by  an  Arbitrary  Cavity -  66 

5-3  Discussion -  72 

5- 4  References  for  Chapter  5 -  73 

CHAPTER  6  ELECTROMAGNETIC  TRANSMISSION  THROUGH  NARROW  SLOTS 

IN  THICK  CONDUCTING  SCREENS -  74 

6- 1  Introduction -  74 

6-2  Formulation  of  the  Problem -  75 

6-3  Specialization  to  Narrow  Slots -  80 

6-4  Numerical  Results -  87 

6-5  Conclusion - 94 

6- 6  References  for  Chapter  6 -  95 

CHAPTER  7  RESONANT  BEHAVIOR  OF  A  SMALL  APERTURE  BACKED  BY  A 

CONDUCTING  BODY -  96 

7- 1  Introduction -  96 

7-2  Formulation  of  the  Problem -  97 

7-3  Aperture  Admittance  Matrix  for  Small  Holes -  104 

7-4  Transmitted  Power -  110 

7-5  An  Example -  113 

7-6  Discussion -  116 

7-7  References  for  Chapter  7 -  118 


iv 


Chapter  I 


BASIC  CONCEPTS 


1-1  INTRODUCTION 

The  coupling  of  electromagnetic  fields  between  isolated  regions 
through  one  or  more  apertures  is  a  widely  encountered  problem  in  elec¬ 
tromagnetics.  Some  examples  are  cavity-to-cavity  coupling,  waveguide- 
to-waveguide  coupling,  waveguide-to-exter ior  space  coupling,  and  so  on. 
The  general  problem  consists  of  two  or  more  regions  of  space  coupled  by 
one  or  more  apertures.  There  can  be  sources  in  one  or  more  regions, 
and  material  bodies  in  one  or  more  regions.  Figure  1-1  shows  a  typical 
problem  of  two  regions  coupled  by  an  aperture.  For  this  example,  it 
is  assumed  that  sources  exist  in  region  a,  and  a  material  body  exists 
in  region  b.  The  boundary  of  the  regions  is  considered  to  be  a  perfect 
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electric  conductor,  but  other  types  of  impermeable  boundaries  may  be 
used.  Each  region  may  be  closed,  as  region  a  is  shown  in  Fig.  1-1,  or 
may  be  open  (extend  to  infinity),  as  region  b  is  shown. 

The  present  chapter  is  intended  to  serve  as  a  short  summary  of 
the  principal  relationships  and  notation  that  we  use  throughout  the 
report.  For  some  readers  this  should  be  only  a  quick  review.  Others, 
less  familiar  with  electromagnetic  theory,  should  refer  to  some  of  the 
well-known  textbooks  [1]  to  [7].  (Bracketed  numbers  denote  references 
which  are  listed  at  the  end  of  each  chapter.) 


1-2  FIELDS  AND  SOURCES 

For  our  purposes,  the  electromagnetic  field  is  expressed  in  terms 
of  two  vectors: 

E  -  the  electric  field 
H  -  the  magnetic  field 

The  source  of  the  field  may  be  of  two  types 

•J  -  the  electric  current 
M  -  the  magnetic  current 

We  use  these  symbols  for  both  volume  densities  and  for  surface  densities. 
It  should  be  clear  from  the  text  or  equations  which  density  is  meant. 

The  medium  jo  which  the  field  exists  is  characterized  by  three  consti¬ 
tutive  parameters. 

£  -  the  permittivity  or  capacitivity 

p  -  the  permeability  or  inductivity 

0  -  the  conductivity 


e  and  u  may  be  complex  to  allow  for  energy  dissipation,  in  which  case  a 

need  not  be  shown  explicitly. 


In  any  region  of  space  for  which  E  and  H  are  differentiable,  the 
field  satisfies  the  Maxwell  curl  equations 


V  *  H  =  jwe  E  +  J 

- V  x  =  jcjy  H  +  M 


(1-1) 


where  V  x  is  the  curl  operator  and  u>  is  the  angular  frequency.  The  diver 
gence  of  (1-1)  yields  the  Maxwell  divergence  equations 


V  •  £E=  q 

V  •  uH= m 

where  V  •  is  the  divergence  operator  and  the  two  scalars  are 


(1-2) 


q  -  the  electric  charge  density 
m  -  the  magnetic  charge  density 

These  are  related  to  the  electric  and  magnetic  currents  by  the  equations 
of  continuity 


V  •  J  =  -  ju>q 

V  •  M  *  -  jam 


(1-3) 


Here  the  symbols  q  and  m  denote  volume  densities,  but  we  use  the  same 
symbols  for  surface  densities.  It  should  be  clear  from  the  text  and 
equations  which  density  is  meant. 

In  any  region  that  is  source-free  (J  =  M  =  0)  ,  linear,  and 
homogeneous  (e  and  M  independent  of  position), the  curl  of  equations 
(1-1)  leads  to  the  vector  Helmholtz  equations 


(1-4) 
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-V  x  V  x  E  +  k  E  =  0 


-VxVxH+kH=0 


where  k  =  w/ey  is  the  wavenumber.  The  rectangular  components  of  IS  and  H 
satisfy  the  scalar  Helmholtz  equation 


V2E,  +  k2E.  =0 
i  i 

(1-5) 

V2Hi  +  k2H±  =  0 

where  i  =  x,  y,  or  z.  In  addition  to  (1-5),  the  vectors  must  satisfy 


V  •  E  =  0  and  V  •  H  =  0. 
1-3  POTENTIALS 


It  is  often  convenient  to  express  the  electromagnetic  fields  in 
terms  of  auxiliary  functions,  called  potentials.  In  a  homogeneous 
region  we  have  the  representation 


E  =  -  jujA  -  V<)> 

H  =  -  V  x  A 
-  y  -  - 

where  the  potentials  are 

A  -  the  magnetic  vector  potential 
<J>  -  the  electric  scalar  potential 
The  potentials  A  and  satisfy  the  Lorentz  gauge 


(1-6) 


V  •  A  =  -  jujye<fr  (1-7) 

which  is  analogous  to  the  equation  of  continuity  (1-3).  Potentials  with 
other  gauges  are  sometimes  useful,  but  we  do  not  consider  them  in  this 


report. 
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Alternatively,  the  field  in  a  homogeneous  region  can  be  represented 
as 

E  =  -  i  y  X  F 

(1-8) 

H  =  -  jcoF  -  Vip 

where  the  potentials  are 

-  the  electric  vector  potential 
>p  -  the  magnetic  scalar  potential 

The  potentials  F  and  ip  also  satisfy  the  Lorentz  gauge 

V  •  F  =  -  jwvieip  (1-9) 

The  potentials  F  and  4'  are  said  to  be  dual  to  the  potentials  A  and  <J>. 

In  general,  the  field  in  a  homogeneous  region  can  be  represented 
as  a  superposition  of  (1-6)  and  (1-8),  or 

E  =  -  i  V  X  F  -  ju)A  -  V(J) 

(1-10) 

H  =  —  V  x  A  -  ju)F  -  V'p 

-  u  -  -  ~  ~ 

In  a  source-free  homogeneous  region,  the  rectangular  components  of  A  and  F 
and  the  scalars  0  and  ip  satisfy  the  Helmholtz  equation 

V2tt  +  k2TT  =  0  (1-11) 

When  one  solution  of  (1-11)  is  taken  for  a  rectangular  component  of  A, 
and  another  solution  of  (1-11)  is  taken  for  a  rectangular  component  of  F, 
the  potentials  are  sometimes  called  Hertzian  potentials. 

One  advantage  of  using  potentials  is  that  the  conditions  V  •  F  =  0 
and  V  •  H  =  0  are  automatically  satisfied.  Another  advantage  is  that  the 
potentials  can  be  simply  related  to  the  sources,  as  we  show  in  the  next 


section. 
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1-4  POTENTIAL  INTEGRALS 

Given  an  electric  current  J_  and  its  associated  electric  change  q 
in  an  infinite  homogeneous  medium  characterized  by  p  and  e,  the  potential 
integrals  are 


(1-12) 


(1-13) 


In  (1-12)  and  (1-13),  i:  is  the  radius  vector  to  the  field  point  and  r_ ' 
is  the  radius  vector  to  the  source  point.  The  sources  J  and  q  are  con¬ 
sidered  to  be  volume  densities  in  (1-12) ,  as  evidenced  by  the  triple 
integral  sign  with  dx '  representing  a  volume  element  of  integration. 

If  the  current  and  charge  are  surface  densities,  it  is  merely  necessary 
to  change  to  a  surface  integration,  and  similarly  for  line  densities. 

The  electromagnetic  field  associated  with  A  and  (J)  is  given  by  (1-6)  . 

The  charge  q  is  related  to  the  current  J  by  the  equation  of  continuity  (1-3). 

For  the  dual  case  of  a  magnetic  current  M  and  its  associated  mag¬ 
netic  charge  m  in  an  infinite  homogeneous  medium  characterized  by  £  and  p, 
the  potential  integrals  are 


F(r)  =  e 


M(r’)  G(r,r')dx' 


m(r')  C(r,r’)dx’ 


(1-14) 
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where  the  Green's  function  G  is  again  given  by  (1-13).  The  electromagnetic 
field  associated  with  £  and  ip  is  given  by  (1-8).  The  charge  m  is  related 
to  the  current  M  by  the  equation  of  continuity  (1-3). 

If  the  problem  of  interest  is  two-dimensional,  that  is,  if  the 
sources  are  infinite  in  extent  in  one  rectangular  direction  and  independent 
of  that  direction,  the  Green's  function  must  be  changed  to 

G(r,  r')  =  H^2)(k|r  -  r' |)  (1-15) 

(2) 

Here  is  the  zero  order  Hankel  function  of  the  second  kind,  equal  to 

Jq  -  jNQ«  where  is  the  zero  order  Bessel  function  of  the  first  kind  and 
Nq  is  the  zero  order  Bessel  function  of  the  second  kind  (or  Neumann  func¬ 
tion).  The  potential  integrals  (1-12)  and  (1-14)  must  now  be  surface 
integrals  in  the  cross  sectional  surface  transverse  to  the  direction  of 
invariance  of  the  current. 

1-5  CONDITIONS  AT  SURFACES 

There  may  be  discontinuities  in  the  electromagnetic  field  at  sheets 
of  current  and  at  surface  discontinuities  of  the  constitutive  parameters. 

In  Fig.  1-2,  let  S  represent  a  surface  between  two  regions,  a  and  b,  and 
let  n  be  the  unit  vector  normal  to  S  pointing  into  region  a.  There  may 
be  an  electric  surface  current  J  and/or  a  magnetic  surface  current  M  on  S. 
(In  figures,  we  use  a  single  headed  arrow  to  denote  electric  current  and 
a  double  headed  arrow  to  denote  magnetic  current.)  The  tangential  com¬ 
ponents  of  the  field  on  S  then  obey  the  conditions 

J  -  n  x  (Ha  -  Hb)  ^ 

M  -  (E3  -  Eb)  x  n 


on  S 


(1-16) 
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Fig.  1-2.  A  boundary  surface  S  between  regions  a  and  b,  possibly 
supporting  surface  currents  J  and  M. 


where  superscript  a  on  a  vector  denotes  that  it  is  evaluated  in  or  is 
peculiar  to  region  a,  and  similarly  for  the  superscript  b.  The  field 
vectors  in  (1-16)  are  evaluated  on  the  surface  S,  as  implied  by  "on  S" 
in  (1-16),  with  the  superscript  denoting  on  which  side  of  S. 

The  normal  components  of  the  fields  on  S  obey  the  conditions 

,  a„a  ~b„b.  ''I 

q  =  n  •  (e  -EE)  I 

V  on  S  (1-17) 

,  a„a  b„b.  1 

m  =  n  •  (u  H  -  U  H  ) 


where  q  and  m  are  the  surface  densities  of  electric  and  magnetic  charges  on 
S,  ea,  u3  are  e,  in  region  a,  and  eb,  y*5  are  e>  y  regi0n  b.  The  q 
and  m  are  related  to  J  and  M  by  the  equation  of  continuity  (1-3),  except 
that  now  V  •  represents  the  surface  divergence  operator. 

If  there  is  no  surface  current  on  S,  then  the  left-hand  sides  of 
(1-16)  are  zero  and 


n  x  (Ha  -  Hb)  =  0 
n  x  (Ea  -  Eb)  =  0 


•  on  S 


(1-18) 
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These  equations  state  that  the  tangential  components  of  ^  and  H  must 
be  continuous  across  S.  There  are  no  surface  densities  of  induced  cur¬ 
rent  at  a  boundary  between  two  media  (perfect  conductors  excepted) . 
Hence,  the  tangential  components  of  _E  and  1!  must  be  continuous  across 
the  boundary  between  two  media  (perfect  conductors  excepted). 

Similarly,  if  there  is  no  surface  charge  on  S,  the  left-hand 
sides  of  (1-17)  are  zero  and 


a_a  b„b 

n  •  e  E  =  n  •  £  E 


on  S 


(1-19) 


a„a  b„b 

n  •  y  H  =  n  *  y  H 


These  equations  state  that  the  normal  components  of  £E  and  pH  must  be 
continuous  across  S.  In  the  time-harmonic  case,  there  are  no  surface 
densities  of  induced  charge  at  a  boundary  between  two  media  (perfect 
conductors  excepted).  Hence,  the  normal  components  of  and  yH  must 
be  continuous  across  the  boundary  between  two  media  (perfect  conductors 
excepted) . 

In  the  special  case  of  a  perfect  electric  conductor,  no  field 
exists  internal  to  the  perfect  conductor,  and  a  surface  current  J  exists 
on  the  surface.  If  region  b  is  a  perfect  conductor,  boundary  conditons 
(1-16)  become 


J  =  n  >:  Ha 
0  =  n  x  Ea 


on  S 


(1-20) 


Hence,  the  tangential  component  of  is  zero  at  the  surface  of  a  perfect 
electric  conductor,  and  the  tangential  component  of  H  is  equal  to  .J  ro¬ 
tated  90°.  Also,  (1-17)  becomes  the  boundary  conditions 


10 


q  =  n 


0  =  n 


a„a 
£  E 


a„a 

Vi  H 


] 


J 


on  S 


(1-21) 


Hence,  the  normal  component  of  pH  is  zero  at  the  surface  of  a  perfect 
electric  conductor,  and  the  normal  component  of  £JS  is  equal  to  the  sur¬ 
face  charge  density. 


1-6  THE  EQUIVALENCE  PRINCIPLE 

Many  source  distributions  outside  a  given  region  can  produce  the 
same  field  within  the  region.  Two  sources  producing  the  same  field  in 
a  region  are  said  to  be  equivalent  within  that  region.  When  we  are 
interested  in  the  field  in  a  given  region,  we  do  not  need  to  know  the 
actual  sources.  Equivalent  sources  serve  as  well. 

A  detailed  discussion  of  the  equivalence  principle  can  be  found 
in  [1] .  Basically,  it  involves  dividing  space  into  two  (or  more)  regions 
and  assuming  Maxwellian  fields  in  each  region.  There  is  generally  a 
discontinuity  in  the  field  at  the  surface  S  which  forms  the  common  boundary 
between  the  two  regions.  On  this  boundary  we  assume  surface  currents  .J 
and  M  which  are  related  to  the  fields  by  (1-16).  We  then  have  a  Maxwellian 
field  everywhere  and  the  sources  which  support  it.  If  there  is  a  one-to- 
one  correspondence  between  a  field  and  its  sources,  we  have  then  found 
equivalent  sources  for  the  field  in  a  given  region.  This  one-to-one 
correspondence  is  always  obtained  in  the  case  of  lossy  media.  It  is  not 
always  obtained  in  the  loss-free  case,  as  noted  in  Section  3-3  of  [1]. 

The  equivalence  principle  is  best  illustrated  by  means  of  a  fairly 


general  example.  Consider  Fig.  1-3  and  let  S  represent  a  mathematical 
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(c)  equivalent  to  a  external  to  S 
and  to  b  internal  to  S 


Fig.  1-3.  Illustration  of  the  Equivalence  Principle. 
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cl  dL 

surface  defined  in  space.  In  Fig.  l-3a,  let  d.  ,  M  be  impressed  sources 

3  3 

producing  a  field  li  ,  H  everywhere  in  a  medium  characterized  by  con¬ 
stitutive  parameters  |ja  and  ea.  In  Fig.  l-3b,  let  J[b,  ^b  impressed 
sources  producing  a  field  E*3,  Hb  everywhere  in  a  medium  characterized  by 
constitutive  parameter j  yb  and  eb.  We  now  construct  the  equivalent 
problem  of  Fig.  3c  as  follows.  External  to  S,  we  specify  that  the  sources 

3  3  3  3  3  3 

are  J  ,  M  ,  the  field  is  _E  ,  1}  ,  and  the  medium  is  y  ,  e  .  Internal  to  S, 
we  specify  that  the  field  is  Eb,  Hb  and  the  medium  is  yb,  £b.  To  support 
the  discontinuity  in  fields  on  S,  we  must  satisfy  (1-16),  or 

J  =  n  x  (H3  -  Hb)  ^ 

V  on  S  (1-22) 

M  =  (Ea  -  Eb)  x  n  j 

These  are  called  the  equivalent  currents  on  S.  Assuming  that  the  relation¬ 
ship  between  the  field  and  its  sources  is  one-to-one,  we  know  that  the 
equivalent  sources  in  Fig.  l-3c  must  produce  the  postulated  field. 

Some  important  special  cases  are:  (1)  the  field  in  one  region  is 
assumed  to  be  the  null  field,  (2)  the  medium  in  one  region  is  assumed  to 
be  a  perfect  electric  conductor,  and  (3)  the  medium  in  one  region  is  assumed 
to  be  a  perfect  magnetic  conductor.  These  cases  are  discussed  in  detail 
in  Section  3-5  of  [1],  and  we  use  such  cases  in  the  next  chapter. 


INTEGRAL  EQUATIONS 


For  some  geometrically  simple  problems,  solutions  for  the  field 
can  be  obtained  by  solving  the  differential  equations.  For  more  compli¬ 
cated  geometries,  it  becomes  more  convenient  to  use  integral  equations 
or,  more  generally  to  use  integro-dif f erential  equations  to  obtain 
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solutions.  These  solutions  are  usually  approximate,  and  one  method  for 
obtaining  them  is  the  method  of  moments,  discussed  in  the  next  section. 

A  general  procedure  for  establishing  integral  equations  is  to  use 
the  potential  integrals  to  find  the  fields  in  terms  of  the  sources  (actual 
or  equivalent)  and  then  to  relate  the  fields  to  the  sources  by  boundary 
conditions  or  constitutive  relationships.  The  procedure  is  best  illus¬ 
trated  by  an  example. 

Let  Pig.  1-4  represent  a  dielectric  body  V  in  free  space,  excited 
by  an  impressed  field  E1,  which  can  be  thought  of  as  the  field  which 
exists  due  to  external  sources  when  the  body  is  absent.  Let  the  con¬ 
stitutive  parameters  of  free  space  be  denoted  by  and  those  of  the 

dielectric  body  by  £,  y^.  The  total  field  within  V  induces  a  polariza¬ 
tion  current  J  related  to  the  total  field  ][  by 

J  =  jco  (e  -  eQ)E  in  V  (1-23) 

In  other  words,  J  is  the  excess  displacement  current  over  what  would 


e‘+  e’ 

€o 


Fig.  1-4.  A  dielectric  body  in  free  space  excited 
by  an  impressed  field  E^. 
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exist  there  if  V  contained  free  space,  E  being  unchanged.  The  total 
field  E  is  the  sum  of  the  impressed  field  E1  plus  the  scattered  (or 

g 

secondary)  field  E  produced  by  .J,  that  is 

E  =  E1  +  ES(J)  (1-24) 

g 

We  can  calculate  IS  (J)  by  the  potential  integrals,  section  1-4,  but 

choose  not  to  write  down  the  explicit  formulas.  For  now,  it  is  suf- 

s  s 

ficient  to  note  that  E  is  linearly  related  to  J,  that  is,  IS  (J)  is  a 

linear  operator.  It  involves  an  integral  and  some  derivatives  and 
should  properly  be  called  an  integro-dif f erential  operator.  For  bre¬ 
vity,  we  shall  call  it  simply  an  integral  operator. 

We  now  use  the  constitutive  relationship  (1-23)  in  (1-24)  to 

obtain 

s  )  -  -1  +  i"  v  O-25) 

J  o' 

This  can  be  rearranged  into  the  form 

L(J)  =  E1  (1-26) 

in  V,  where 

'  WT  1  )  -  I8®  f1"27) 

J  o 

is  a  linear  operation.  In  (1-26),  the  impressed  (or  incident)  field  E1 
is  known,  L  is  a  known  integral  operator,  and  J  is  the  unknown  to  be 


found. 
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1-8  METHOD  OF  MOMENTS 

The  method  of  moments  is  a  general  procedure  for  solving  operator 
equations  of  the  form  (1-26).  Its  mathematical  foundations  lie  in  the 
theory  of  projections  in  linear  inner  product  spaces.  A  detailed  but 
elementary  exposition  of  the  procedure  can  be  found  in  reference  [8], 

The  following  is  a  short  outline  of  the  method. 

Suppose  we  have  an  operator  equation  of  the  form 

L(f)  =  g  (1-28) 

where  L  is  a  linear  operator,  g  is  a  known  function,  and  f  is  an  unknown 
function  to  be  determined.  We  must  define  an  inner  product  <f,  g>  (or 
symmetric  product  if  f  and  g  are  complex)  such  that 

(a)  <f,  g>  =  <g,  f> 

(b)  <af1  +  8f2,  g>  =  a<flf  g>  +  0<f2,g>  (1-29) 

(c)  <f,  f*>  >  0  if  f  i*  0 

(d)  <f,  f  >  =  0  if  and  only  if  f  =  0 

Here  a  and  8  are  scalars,  and  *  denotes  complex  conjugate.  In  a  complex 

* 

Hilbert  space  the  usual  inner  product  is  <f,  g  >,  and  the  conjugate  oper¬ 
ation  is  not  shown  explicitly  .  In  electromagnetic  theory  it  is  more  con¬ 
venient  to  use  the  symmetric  product  <f,  g>  and  show  the  conjugation 
explicitly  when  needed. 

We  next  choose  a  set  of  linearly  independent  expansion  functions 
or  basis  functions  { f ^ ,  f^,...}  in  the  domain  of  L  and  represent  f  as  the 


linear  combination 
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f 


*  l 


a  £ 
n  n 


(1-30) 


where  the  are  scalars  to  be  determined.  For  computational  purposes 
the  representation  is  usually  approximate,  as  shown.  Substituting 
(1-30)  into  (1-28)  and  using  the  linearity  of  L,  we  obtain 


(1-31) 


We  next  define  a  set  of  linearly  independent  testing  functions  or  weight¬ 
ing  functions  {w^,  w^,...)  in  the  range  of  L  and  take  the  inner  product 

(or  symmetric  product)  of  (1-31)  with  each  w  .  The  inner  product  is 

m 

linear,  and  the  result  is 


l  a  <w  ,  Lf  >  =  <w  ,  g>  (1-32) 

n  m  n  m 

n 

m  =  1,2,3,...  .  This  is  a  set  of  linear  equations  which  can  be  written 
in  matrix  form  as 

[A]a  =  g  (1-33) 

where 


m  - 

<w1,  Lf1> 

<W1’ 

Lf 2>  •  •  • 

<w2,  Lf1> 

<w2. 

Lf 2>  .  .  . 

•  *  * 

•  •••«« 

• 

* 

°1 

g  = 

1’ 

g> 

a2 

<w2, 

g> 

• 

(1-34) 


(1-35) 


If  the  number  of  expansion  functions  equals  the  number  of  testing  functions 
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[£]  will  be  a  square  matrix.  The  a  and  g  are  column  matrices,  also  called 
column  vectors.  If  [£]  is  nonsingular,  its  inverse  [£]  ^  exists,  and  the 
a  is  then  given  by 

a  =  m-1g  (1-36) 

The  solution  for  f  is  then  given  by  (1-30)  where  a  are  the  components  of 
a  obtained  from  (1-36) . 

By  defining  a  row  matrix  of  the  expansion  functions  as  (where 
the  tilde  is  used  to  denote  the  transpose  of  a  matrix  In  general) 

f  =  [f-L  f2  f3  •  •  •]  (1-37) 

we  can  write  the  solution  concisely  in  the  form 

f  ~  ?  a  =  fU]-1!  (1-38) 

Finally,  instead  of  the  function  f  itself,  we  are  often  interested  in 
some  functional  (number)  p  which  depends  linearly  on  f.  This  can  be 
expressed  in  the  general  form 

p  =  <h,  f>  (1-39) 

where  h  is  a  known  function  (determined  as  a  part  of  formulating  the 
problem).  Substituting  (1-38)  into  (1-39),  we  have 

P  *  h[ £] _1g  (1-40) 

where  h  is  the  row  vector 

h  =  [<f1>  h>  <f2,  h>  ...]  (1-41) 

In  general,  we  call  h  a  measurement  vector.  (It  is  characteristic  of 
the  particular  measurement  p  we  wish  to  perform.)  Every  linear  measurement 
(functional)  can  be  expressed  as  a  matrix  contraction  of  the  form  (1-40). 
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Chapter  II 

A  GENERALIZED  NETWORK  FORMULATION  FOR  APERTURE  PROBLEMS 

2-1  INTRODUCTION 

This  chapter  considers  a  general  formulation  for  a  two-region  aper¬ 
ture  problem.  First  an  operator  equation  is  obtained  in  terms  of  an 
unknown  equivalent  magnetic  current,  and  this  is  then  reduced  to  a  matrix 
equation  via  the  method  of  moments.  The  only  coupling  is  through  the 
aperture,  whose  characteristics  are  expressed  by  aperture  admittance 
matrices,  one  for  each  region.  These  admittance  matrices  depend  only  on 
the  region  being  considered,  being  independent  of  the  other  region.  The 
aperture  coupling  is  then  expressible  as  the  sum  of  the  two  independent 
aperture  admittance  matrices,  with  source  terms  related  to  the  incident 
magnetic  field.  This  result  can  be  interpreted  in  terms  of  generalized 
networks  as  two  N-port  networks  connected  in  parallel  with  current 
sources.  The  resultant  solution  is  equivalent  to  an  N-term  variational 
solution. 

Since  the  problem  is  divided  into  two  mutually  exclusive  parts,  one 
can  separately  solve  a  few  canonical  problems,  such  as  apertures  in  con¬ 
ducting  screens,  in  waveguides,  and  in  cavities,  and  then  combine  them  in 
various  permutations.  Computer  programs  can  be  developed  to  calculate 
the  aperture  admittance  matrices  for  classes  of  canonical  problems, 
such  as  apertures  of  arbitrary  shape  in  conducting  planes,  in  square  wave¬ 
guides,  and  in  rectangular  cavities.  Such  programs  can  then  serve  as 
broad  and  versatile  tools  for  designing  electromagnetic  networks  with  aper¬ 


ture  coupling. 


2-2  GENERAL  FORMULATION 

Figure  2-1  represents  the  general  problem  of  aperture  coupling  between 
two  regions,  called  region  a  and  region  b.  In  region  a  there  are  impressed 
sources  J_  ,  M*,  and  region  b  is  assumed  source  free.  The  more  general  case 
of  sources  in  both  region  a  and  region  b  can  be  treated  as  the  superposition 
of  two  problems,  one  with  sources  in  region  a  only,  plus  one  with  sources  in 
region  b  only.  Each  region  of  Fig.  2-1  is  shown  to  be  bounded  by  an  elec¬ 
tric  conductor,  although  other  types  of  electromagnetic  isolation  may  be 
used.  Region  a  is  shown  closed  and  region  b  is  shown  open  to  infinity, 
although  each  region  may  be  open  or  closed.  The  equivalence  principle, 
discussed  in  Sec.  1-6,  is  used  to  divide  the  problem  into  two  equivalent 


Fig.  2-1.  The  general  problem  of  two  regions 
coupled  by  an  aperture. 
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problems  as  shown  in  Fig.  2-2.  In  region  a,  the  field  is  produced  by 
the  sources  J*,  M*,  plus  the  equivalent  magnetic  current 

M  =  n<  E  (2-1) 

over  the  aperture  region,  with  the  aperture  covered  by  an  electric  con¬ 
ductor.  In  region  b,  the  field  is  produced  by  the  equivalent  magnetic 
current  -M  over  the  aperture  region,  with  the  aperture  covered  by  an 
electric  conductor.  The  fact  that  the  equivalent  current  in  region  b  is 
the  negative  of  that  in  region  a  ensures  that  the  tangential  component  of 
electric  field  is  continuous  across  the  aperture.  The  remaining  boundary 
condition  to  be  applied  is  continuity  of  the  tangential  component  of  mag¬ 
netic  field  across  the  aperture. 

The  tangential  component  of  magnetic  field  in  region  a  over  the 

aperture,  denoted  ,  is  the  sum  of  that  due  to  the  impressed  sources, 

i  3 

denoted  H  ,  plus  that  due  to  the  equivalent  source  M,  denoted  H  (M) , 

— t  -  — t  — 

that  is 

H^  =  Hj  +  H^(M)  v 

1.  3 

Note  that  H  and  H  (M)  are  both  computed  with  a  conductor  covering  the 
— t  — t  — 

aperture.  A  similar  equation  holds  for  region  b,  except  that  the  equi¬ 
valent  source  -M  is  the  only  source.  Hence,  the  tangential  component 
of  magnetic  field  in  region  b  over  the  aperture  is 

H*  =  H?(-M)  =  -  H?(M)  (2-3) 

t  t  —  t  — 

where  H^(M)  is  computed  with  a  conductor  covering  the  aperture.  The  last 
equality  in  (2-3)  is  a  consequence  of  the  linearity  of  the  operator. 
The  true  solution  is  obtained  when  nf  of  (2-2)  equals  Hb  of  (2-3),  or 


CONDUCTOR 


(a)  EQUIVALENCE  FOR  REGION  a. 


(b)  EQUIVALENCE  FOR  REGION  b. 

Fig.  2-2.  The  original  problem  divided  into  two 
equivalent  problems. 
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-  h!(M)  -  H^(M)  =  H*  (2-4) 

— t  —  “t  ~  ~t 

This  is  the  basic  operator  equation  for  determining  the  equivalent 
magnetic  current  M. 

If  (2-4)  were  satisfied  exactly,  we  would  have  the  true  solution. 

We  use  the  method  of  moments  to  obtain  an  approximate  solution. 
Define  a  set  of  expansion  functions  (M  ,  n=l,2, . . . ,N},  and  let 


<A,B>  = 


A  •  B  ds 


(2-7) 


apert. 


and  a  set  of  testing  functions  {W^  . . . ,N),  which  may  or  may  not  be 
equal  to  the  expansion  functions.  We  take  the  symmetric  product  of  (2-6) 
with  each  testing  function  W  ,  and  use  the  linearity  of  the  symmetric 


product  to  obtain  the  set  of  equations 


m»l,2, . . . ,N.  Solution  of  this  set  of  linear  equations  determines  the 

coefficients  V  and  the  magnetic  current  M  according  to  (2-5).  Once  M 
n  —  — 

is  known,  the  fields  and  field-related  parameters  may  be  computed  by 


standard  methods. 
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The  above  solution  can  be  put  into  matrix  notation  as  follows: 
Define  an  aperture  admittance  matrix  for  region  a  as 


[Ya]  =  [<-W  ,Ha(M  )>1 

m  t  NxN 


and  an  aperture  admittance  matrix  for  region  b  as 


[Yb]  =  [<-W  ,Hb(M  )>} 

m  t  NxN 


(2-9) 


(2-10) 


The  minus  signs  are  placed  in  (2-9)  and  (2-10)  on  the  basis  of  power  con¬ 
siderations.  Define  a  source  vector 


^  ‘  'V^Nxl 


(2-11) 


and  a  coefficient  vector 


V  ’  'V»xl 


(2-12) 


Now  the  matrix  equation  equivalent  to  equations  (2-8)  is 


[Ya  +  Yb]V  =  I1 


(2-13) 


This  can  be  interpreted  in  terms  of  generalized  networks  as  two  networks 


Fig.  2-3.  The  generalized  network  interpretation 
of  equation  (13). 
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[Ya]  and  [Y^]  in  parallel  with  the  current  source  I*,  as  shown  in  Fig.  2-3. 
The  resultant  voltage  vector 

V  =  [Ya  +  Y15]'1!1  (2-14) 

is  then  the  vector  of  coefficients  which  gives  M  according  to  (2-5) . 

g 

It  is  important  to  note  that  computation  of  [Y  ]  involves  only 
region  a,  and  computation  of  fY^]  involves  only  region  b.  Hence,  we 
have  divided  the  problem  into  two  parts,  each  of  which  may  be  formulated 
independently.  Once  [Ya]  is  computed  for  one  region,  it  may  be  combined 

g 

with  [Y  ]  for  any  other  region,  making  it  useful  for  a  wide  range  of 
problems.  For  example,  the  same  aperture  admittance  matrix  for  radi¬ 
ation  into  half-space  would  be  useful  for  plane  wave  excitation  of  the 
aperture,  waveguide  excitation,  and  cavity  excitation. 

2-3  LINEAR  MEASUREMENT 

A  linear  measurement  is  defined  as  a  number  which  depends  linearly 
on  the  source.  Examples  of  linear  measurements  are  components  of  the  field 
at  a  point,  voltage  along  a  given  contour,  and  current  crossing  a  given 
surface.  Measurements  made  in  region  b  will  depend  linearly  only  on  the 
equivalent  current  -M.  Measurements  made  in  region  a  will  depend  linearly 
on  the  impressed  sources  J*',  M^,  as  well  as  on  the  equivalent  current  M. 

We  now  illustrate  these  concepts  with  a  particular  example. 

Consider  the  measurement  (computation)  of  a  component  ii^  of  mag¬ 
netic  field  at  a  point  r^  in  region  b.  It  is  known  that  this  component 
can  be  obtained  by  placing  a  magnetic  dipole  at  r^,  and  applying  the 
reciprocity  theorem  to  its  field  and  to  the  original  field.  Section  3-8 
of  [1].  The  original  field  in  region  b  is  given  by  the  solution  to 
Fig.  2-2b.  The  problem  involving  the  magnetic  dipole,  called  the  adjoint 
problem,  is  shown  in  Fig.  2-4.  Application  of  the  reciprocity  theorem 
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Fig.  2-4.  The  adjoint  problem  for 

determining  H  at  r  . 

m  — m 


to  these  two  cases  yields 


H  K1  =  - 
m  m 


M  •  Hmds 


(2-15) 


apert. 


Here  H  is  the  magnetic  field  from  K£ ^  in  the  presence  of  a  complete 

conductor,  and  H  is  the  component  in  the  direction  of  K£  of  the  mag- 
m  -m 

netic  field  at  r^  due  to  -M  in  the  presence  of  a  complete  conductor. 

To  evaluate  (2-15),  substitute  for  M  from  (2-5)  and  obtain 


H  KS2  =  Y  V  <-M  ,  Hm> 
m  m  L  n  n 


(2-16) 


This  can  be  written  in  matrix  form  as 


H  KK  =  I  V 
m  m 


where  1  is  the  transpose  of  a  measurement  vector 


(2-17) 
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I®  =  [<-M  ,  Hm> ] 

n  Nxl 


(2-18) 


Note  that  the  elements  of  I™  are  similar  in  form  to  those  of  I*  given 

by  (2-11),  except  that  -M  replaces  W  .  The  minus  sign  difference 

— n  ~ m 

reflects  the  fact  that  the  equivalent  source  in  region  b  is  -M,  in  con¬ 
trast  to  that  in  region  a  which  is  +M.  Now  substitute  (2-14)  into  (2-17) 
to  obtain 


HKl  =  lm[Ya  +  Yb]-1  I1 
mm 


(2-19) 


If  the  magnetic  dipole  is  of  unit  moment,  then  (2-19)  gives  H  at  r 

m  — m 

directly. 

Every  linear  measurement  in  region  b  will  be  of  the  form  (2-19). 

For  example,  if  a  component  of  E  at  were  desired,  we  would  place  an 
electric  dipole  at  r^  and  apply  reciprocity.  In  general,  a  linear  measure¬ 
ment  involves  applying  reciprocity  to  the  original  problem  and  to  an  ad¬ 
joint  problem.  A  determination  of  the  sources  of  the  adjoint  problem  is  a 
part  of  the  formulation  of  the  problem. 

If  a  linear  measurement  is  made  in  region  a,  it  will  involve  a  con¬ 
tribution  from  the  impressed  sources  M*  added  to  that  from  the  equivalent 
source  M.  For  example,  instead  of  (2-19)  we  would  have 


HKl  '  H^fc  +  im[Ya  +  Yb]-1T 
m  m  m  m 


(2-20) 


where  H*  is  the  magnetic  field  from  J^,  M*  in  the  presence  of  a  complete 
m  —  — 

conductor.  Also,  in  region  a  we  would  define  the  measurement  vector  to  be 


T  -  [<M  ,  Hm>]M 

n  Nxl 


(2-21) 


instead  of  (2-18),  because  the  equivalent  sources  are  +M  in  region  a  in 


28 


contrast  to  -M  in  region  b.  Note  that  it  is  the  difference  field  H-H* 
in  region  a  (due  to  M)  that  is  directly  analogous  to  the  transmitted 
field  H  in  region  b  (due  to  -M)  . 


2-4  TRANSMITTED  POWER 

A  quadratic  measurement  is  one  which  depends  quadratically  on  the 
sources.  Examples  of  quadratic  measurements  are  components  of  the 
Poynting  vector  at  a  point,  power  crossing  a  given  surface,  and  energy 
within  a  given  region.  A  particular  quadratic  measurement  of  considerable 
interest  is  the  power  transmitted  through  the  aperture,  which  we  now 
consider. 

The  complex  power  transmitted  through  the  aperture  is  basically 

P  =  E  x  H*  •  n  ds  (2-22) 

t  jj  ~  ~  ~ 

apert. 

where  the  asterisk  denotes  complex  conjugate.  Substituting  from  (2-1), 
we  have 

P  -  [[  M  •  H*ds  (2-23) 

t  JJ  -  - 

apert . 

This  involves  only  the  tangential  component  of  H,  which  in  region  b  we 
denoted  by  H^(-M).  For  M  we  use  the  linear  combination  (2-5)  and  ob¬ 
tain 

H*?(-M)  =  -  7  V_Hb(M  )  (2-24) 

t  —  n  t  ~"n 

n 

Substituting  this  for  H  and  (2-5)  for  M  into  (2-23),  we  obtain 


P_  =  ~  I  l  v  v*  M  •  H^M  )ds  (2-25) 

t  4-i'mn  — m  — t  — n 

m  n  •* 


apert . 


If  M  are  real,  the  conjugate  operations  can  be  taken  outside  the 
— m 
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if  u  -  u  (Galerkin’s  method),  then  the  negative 
integrals.  Moreover,  If  M„  \  <-baierK 

i  in  (?  251  are  Y*3*  as  defined  by  (2-10),  and 
of  the  integrals  in  (2-25)  are  ^  a 

*  b* 

P  =  y  y  v  v  y° 

t  L  L  m  n  mn 

u  m  n 


(2-26) 


rhis  can  be  written  in  matrix  form  as 


p  =  v[vDl  v 


(2-27) 


Sote  that  this  is  the  usual  formula  for  power  into  network  [Y  ]  of 
Fig.  2-3. 


2-5  DISCUSSION 

The  basic  focal.-, tier,  for  two-region  aperture  problems  is  given 
in  this  chapter.  Only  a  single  aperture  is  considered  explicitly,  hut 
the  extension  to  multiple  apertures  is  straightforward .  In  Chapter  Ill. 
we  apply  the  formulation  to  the  problem  of  apertures  In  plane  conducting 
screens.  In  Chapter  IV,  we  apply  it  to  waveguide-fed  apertures.  In 
Chapter  V,  the  formulation  is  specialised  to  cavity-backed  apertures. 
Chapter  VI  considers  the  case  of  a  narrow  slot  in  a  thick  conducting 
plane.  Chapter  VII  considers  the  case  of  an  aperture  in  a  plane  screen 
hacked  by  a  conducting  body.  Chapter  VIII  considers  the  ecuivalent 
circuit  for  coupling  through  an  aperture  to  a  long  wire.  There  arc, 
of  course,  infinitely  many  geometries  that  can  be  considered,  and  this 
report  treats  only  a  representative  cross  section  of  these  possibilities. 
The  basic  principles  remain  the  same  regardless  of  the  type  of  problem 
being  considered. 
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Chapter  III 


APERTURES  IN  PLANE  CONDUCTORS 

3-1  GENERAL  THEORY 

Consider  a  conducting  plane  covering  the  z=0  plane  except  for  an 
aperture,  as  shown  in  Fig.  3-1.  The  two  regions  z>0  and  z<  0  are 
identical-half  spaces,  and  hence  their  aperture  admittance  matrices 
are  the  same.  Therefore,  we  let 

[Ya  +  Yb]  =  2[Yhs]  (3-1) 

tis 

where  [Y  ]  denotes  the  aperture  admittance  for  the  aperture  opening 
into  half  space,  say  z  >  0.  When  the  aperture  is  covered  by  a  conductor. 


Fig.  3-1.  Aperture  in  a  plane  conductor. 
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the  z  =  0  plane  is  a  complete  conducting  plane,  and  image  theory  applies. 

The  magnetic  current  expansion  functions  are  on  the  surface  of  the  z=0 

plane.  Their  images  are  equal  to  them  and  are  also  on  the  z  =  0  plane, 

hs 

according  to  section  3-6  of  [1].  The  result  is  that  [Y  ]  is  the  admit¬ 
tance  matrix  obtained  using  expansion  functions  2M  radiating  into  free 

— n 

space  everywhere.  This  problem  is  dual  to  that  for  the  impedance  matrix 

of  a  plane  conductor,  a  problem  considered  recently  in  the  literature  [2]. 

The  original  excitation  of  the  aperture  is  by  the  impressed 

sources  J/L,  M1  in  the  region  z  <  0.  The  impressed  field  H*  used  in  the 

~t 

operator  equation  (2-4)  is  the  tangential  magnetic  field  due  to  J^,  M"*" 
with  the  aperture  covered  by  a  conductor  (Fig.  2-2a) .  In  this  case  the 
z  =  0  plane  is  a  complete  conductor,  and  image  theory  again  applies.  The 
result  is  that  the  tangential  component  of  H  over  the  z  =  0  plane  when  it 
is  covered  by  a  conductor  is  just  twice  what  it  is  for  the  same  sources 
in  free  space.  Hence, 


2H 


io 


(3-2) 


,io 


where  H  '  is  the  tangential  component  of  the  magnetic  field  over  the 
aperture  due  to  the  sources  J*,  M1  in  free  space.  The  components  of 
the  excitation  vector  I  defined  by  (2-11)  are  now 


1=2 

m 


W  •  Hlods 


(3-3) 


apert . 

where  W  is  the  mth  testing  function. 

A  case  of  special  interest  is  that  of  plane  wave  excitation. 

A  unit  plane  wave  is  given  by 

.  -jk,  •  r 

H  =  u^  e 


(3-4) 
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where  is  a  unit  vector  specifying  the  direction  of  H10,  k_  is  the 
propagation  vector  of  magnitude  2tt/\  and  pointing  in  the  direction  of 
propagation,  and  r_  is  the  radius  vector  to  an  arbitrary  field  point. 
These  vectors  are  shown  in  Fig.  3-1.  The  components  (3-3)  of  the 
plane-wave  excitation  vector  are  then 


P 


i 


m 


apert . 


W 

—m 


(3-5) 


The  symbol  P  has  been  used  for  this  particular  vector  to  distinguish 

it  from  the  more  general  excitation  vector  (3-3). 

Similar  simplifications  apply  to  the  adjoint  (measurement) 

problem.  For  the  evaluation  of  a  component  of  magnetic  field  at  a 

point  r  ,  a  magnetic  dipole  is  placed  at  the  measurement  point  r  . 

~ m  — m  — m 

This  radiates  in  the  presence  of  a  complete  conductor  over  the  z  =  0 
plane,  and  hence,  analogous  to  (3-2),  we  have 

H™  =  2jf°  (3-6) 

Here  H^  denotes  the  tangential  component  of  11  over  the  aperture  from 

when  the  z  =  0  plane  is  covered  by  a  conductor,  and  H™°  denotes  that 

from  K2  when  it  radiates  into  free  space.  The  components  of  the 
~*rn 

measurement  vector  I  defined  by  (2-18)  are  now 

Im  =  -  2  ff  M  •  Hm°ds  (3-7) 

n  JJ  - n  -t 

apert . 

where  is  the  nth  expansion  function. 

A  case  of  special  interest  is  that  of  far-field  measurement. 

This  is  obtained  by  a  procedure  dual  to  that  used  for  radiation  and 
scattering  from  conducting  wires  [3].  To  obtain  a  component  of  H 
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on  the  radiation  sphere,  we  take  a  source  K£  perpendicular  to  r^ 

and  let  r  ■>  At  the  same  time  we  adjust  K£  so  that  it  produces 
m  m  r 

a  unit  plane  wave  in  the  vicinity  of  the  origin.  The  required 
dipole  moment  is  given  by 


,  •  -  j  kr 

1  =  -JGJC  J  m 

K£  47Tr  6 

m  m 


(3-8) 


and  the  plane  wave  field  it  produces  in  the  vicinity  of  the  origin  is 


-jk  •  r 
„mo  — m  — 

H  =  u  e 
—  m 


(3-9) 


Here  u  is  a  unit  vector  in  the  direction  of  H  ,  k  is  the  propaga- 
-m  —  — m 

tion  vector,  and  £  is  the  radius  vector  to  an  arbitrary  field  point. 
Again  these  vectors  are  shown  in  Fig.  3-1.  The  components  (3-7)  of  the 
far-field  measurement  vector  are  then 


P  =  -  2 
n 


-jk  •  £ 

M  •  u  e  ^ 01  ds 

— n  — m 


(3-10) 

11  IU 

apert . 

The  symbol  P'“  is  used  for  this  particular  measurement  vector  to 

distinguish  it  from  the  more  general  measurement  vector  (3-7).  The 

far-zone  magnetic  field  is  now  given  by  (2-19)  with  K£  given  by  (3-8), 

m 

-Kn  -H  <3  K 

I  =  P  ,  I  =  P  ,  and  [Y  +  Y  ]  given  by  (3-1) .  Hence 


— 1  Ur 

H  =  ^e  mpm[YhSr1p1 
m  otr 

m 


(3-11) 


The  usual  two  radiation  components  Hg  and  are  obtained  by  orienting 
K?  in  the  0  and  j)  directions,  respectively. 
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3-2  TRANSMISSION  PARAMETERS 

A  parameter  sometimes  used  to  express  the  transmission  character¬ 
istics  of  an  aperture  is  the  transmission  cross  section  t.  It  is  de¬ 
fined  as  that  area  for  which  the  incident  wave  contains  sufficient  power 

to  produce  the  radiation  field  H  by  omnidirectional  radiation  over  half 

m 

space.  For  unit  incident  magnetic  field,  this  is 


T  = 

Substituting  from  (3-11), 


2nr  H 
m1  m 


we  obtain 


X 


2  2 
a  £ 

32it  |r 


[Yhs]_1  p1! 


(3-12) 


(3-13) 


Note  that  x  depends  upon  the  polarization  and  direction  of  the  incident 
wave  (via  P^) ,  and  upon  the  polarization  measured  and  direction  to  the 
measurement  point  (via  P  ) . 

Another  parameter  used  to  express  the  transmission  characteristics 
of  an  aperture  is  the  transmission  coefficient  T,  defined  as 


T 


trans. 


inc. 


(3-14 ) 


where  P  is  the  time-average  power  transmitted  by  the  aperture, 

trdns • 

and  P  is  the  free  space  power  incident  on  the  aperture.  The  inci- 
me  • 

dent  power  is 


P 


inc. 


*  n  s  cos  9 . 

inc . 


(3-15) 


where  r)  »  /y/ e  is  the  intrinsic  impedance  of  free  space,  S  is  the 

aperture  area,  and  9.  is  the  angle  between  k.  and  n.  The  trans- 

inc  •  i 

mitted  power  is 


3  b 


P  =  Re(P  ) 

trans  t 


(3-16) 


where  Re(Pt)  denotes  the  real  part  of  P^,  given  by  (2-27)  so  that 


T=n~S  core—  Re(VtYhs]*  v*)  (3-17) 

inc . 

Note  that  T  depends  on  both  the  direction  of  incident.  ■  and  on  the 
polarization  of  the  incident  wave. 

Finally,  because  of  symmetry  about  the  z  =  0  plane,  the  difference 
field  H-H1  which  exists  in  the  region  z<0  is  simply  related  to  the 
transmitted  field  which  exists  in  the  region  z>0.  The  difference  field 
in  the  region  z  <  0  is  produced  by  an  equivalent  current  M  on  a  plane 
conductor  over  the  z  =  0  plane.  By  image  theory,  it  is  also  the  field 
produced  in  the  region  z<0  by  the  source  2M  in  free  space.  Analogously, 
the  transmitted  field  in  the  region  z>0  is  produced  by  the  source  -2M 
in  free  space.  Hence,  the  difference  field  in  the  region  z<  0  and  the 
negative  of  the  transmitted  field  in  the  region  z>0  are  both  pro¬ 
duced  by  the  same  magnetic  current  2M  radiating  in  free  space. 


3-3  ADMITTANCE  MATRIX 

If  an  admittance  matrix  [Y]  is  defined  by  [Y]  =  [Y3  +  Y*3]  ,  then, 
according  to  (3-1)  and  (2-10),  the  ij -th  element  of  [Y]  is  given  by 


Yij  =  (y3  +  Y\j  =  '  4<Wi’  H(^j)> 


(3-18) 


where  II (Mj )  is  the  magnetic  field  produced  by  radiating  in  free 
space.  In  view  of  Sections  1-3  and  1-4,  the  magnetic  field  H^(M^) 
can  be  expressed  in  terms  of  an  electric  vector  potential  Fj  and  mag¬ 
netic  scalar  potential  ipj  as  [4] 
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where 


H(M.)  =  -  jojFj  -  Vi|>j 


F  =  ~ 


e_ 

4tt 


-jk!r-r' 


1 1 

apert . 


M.  - 1 - ; 

-3  |r-r' 


ds' 


= 


1 

4ttp 


e-jk|r-r’ 

”3  k-r'T 


ds ' 


apert. 


(3-19) 


(3-20) 


(3-21) 


(3-22) 


where  r  and  _r'  are  respectively  the  vectors  to  the  field  and  source 
points  in  the  aperture.  Substituting  (2-7)  and  (3-19)  into  (3-18),  we 
obtain 


rr 


*«-*  i 


W,  •  (jwF .  +  Vv,)ds  .  (3-23) 

-i  J  ~J  ~  J 

apert . 

If  the  component  of  normal  to  the  rim  of  the  aperture  vanishes 
on  the  rim  of  the  aperture,  then 

V  •  (0>.  W.)  ds  =  0 

J  J  J  1 

apert. 

The  above  equation  can  be  rewritten  as 


J 


W.  •  VO),  ds  + 


<t>  V  •  W.  ds  =  0 
J  ~ 


apert . 

Hence,  (3-23)  becomes 


apert . 


Y«  • 4j" 


(F  •  W.  +  il’.m.)  ds 
-j  -i  "j  i 


where 


apert . 


V  •  W 

«.  =  - ; - 

i  -jw 


(3-24) 


(3-25) 


( 3-26) 


| 
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We  must  now  consider  a  specific  problem  in  order  to  choose  appropriate 
expansion  and  testing  functions. 

3-4  THE  RECTANGULAR  APERTURE 

The  geometry  and  coordinate  system  for  the  rectangular  aperture 
in  a  conducting  plane  is  shown  in  Fig.  3-2.  For  this  problem,  we 
choose  the  set  of  testing  functions  W.  equal  to  the  set  of  expansion 


Fig.  3-2.  Rectangular  aperture  in  a  conducting 
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where  L  and  L  are  integers  is  divided  into  rectangular  subareas 
x  y 

of  length  Ax  in  x  and  Ay  in  y.  The  set  of  expansion  functions 
is  split  into  a  set  of  x  directed  magnetic  currents  and  a  set 
of  y  directed  magnetic  currents  defined  by 


^+(q-l)(Lx-l)  TPW  Pq(j,K 


p  =  1,2, . . .L  -1 
x 

q  =  1,2,.. .  L 


(3-27) 


■  i  ^  Pp« 


p  "■  1 , 2  ,  .  .  .  L 

x 

q  =  1,2, ...L-l 


where  x  and  £  are  unit  vectors.  T*(x)  and  T^(y)  are  triangle 
functions  defined  by 


(3-28) 


T*(x) 

P 


(p-1)  Ax  <  x  pAx 
pAx  £  x  £  (p+l)Ax 
]  x  -  pAx|  >  Ax 


(3-29) 


T^(y)  = 


(q-l)Ay  1  y  1  qAy 
qAy  <  y  <_  (q+l)Ay 


[y  -  qAy |  >  Ay 


(3-30) 


X  V 

and  Pp(x)  and  P^(y)  are  pulse  functions  defined  by 


P*(x)  = 


1  (p-1) Ax  <  x  <  pAx 


0  all  other  x 


(3-31) 


P^(y) 


(q-l)Ay  <  y  <  qAy 


all  other  y 


(3-32) 
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xy  x  y 

The  magnetic  charge  sheets,  say  itk  and  m^  associated  with  and  M' 

are  obtained  from  (3-22)  as 


P+(q-l)(Lx-l) 


m 


(pp(x)  -  ppti(x))<(y) 
-jcuAx 


<pp<p>  - 


P+(q-l)Lx  -jwAy 


(3-33) 


(3-34) 


Introduction  of  the  two  types  of  expansion  functions  and  and  the 
two  types  of  testing  functions  and  into  (3-25)  gives  rise  to 
four  Y  submatrices  defined  by 

f u  =  x,y 


YUJ  =  4JU 


JJ 

apert . 


(F^  •  m“  +  i|>V?)ds 

“j  ~i  J  i 


(3-35) 


x,y 


The  mathematical  details  and  approximations  for  numerically  evaluating 
(3-35)  can  be  found  in  a  research  report  [5]. 


3-5  PLANE  WAVE  EXCITATION  AND  MEASUREMENT  VECTORS 

The  plane  wave  excitation  vector  P^  of  (3-5)  and  the  plane  wave 
measurement  vector  P™  of  (3-10)  are  of  the  same  form  except  for  a  minus 
sign.  We  therefore  need  to  evaluate  only  one  of  them,  say  the  measure- 

-►tn 

ment  vector  P  .  We  specialize  it  to  four  principal  plane  patterns  as 


mu 

^  i  '9y 


=  -  2 


J 

apert . 


Uu  «  Ikx  cos  0  ,  , 

•  0_  e  dxdy 


(P”U) 

i  yy 


-  2 


)) 


M  •  £  e 


jkx  cos  0 


dxdy 


(3-36) 


(3-37) 


apert. 
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,mus  =  _  ? 

i  ;<j>x 

* 

. 

* 

-i 

~  jky  cos 

•  (ti  J 

dxdy 

(3-38) 

apert . 

'tj1U)  =  -  2 

• 

' 

m" 

^  iky  cos 
•  x  J 

<P 

dxdy 

(3-39) 

i  xx 

. 

“"I 

apert. 


where  u  is  either  x  or  y.  The  superscript  u  is  necessary  because  has 
been  split  up  into  and  of  (3-27)  and  (3-28).  In  (3-36)  to  (i-39) , 

£,  and  x  are  unit  vectors  in  the  0,  y,  <j>,  and  x  directions  respec¬ 
tively  where,  as  shown  in  Fig.  3-2,  0  is  measured  from  the  positive  x  axis 
in  the  y=0  plane  and  <f>  is  measured  from  the  positive  y  axis  in  the  x=0 
plane.  For  measurement  vectors,  0°  £  0  £  180°,  0°  £  <f>  £  180°.  (P^U)^  is 
for  a  polarized  measurement  in  the  y=0  plane,  (p™U)yy  is  f°r  a  Z 

polarized  measurement  in  the  y=0  plane,  (p™U)(j)x  is  f°r  a  i  polarized 
measurement  in  the  x=0  plane,  and  (P™11)  is  for  an  x  polarized  measure- 

1  XX 

ment  in  the  x=0  plane.  Because  our  set  of  testing  functions  W4  is  the 

same  as  the  set  of  expansion  functions  ,  the  plane  wave  excitation 

vector  P*  of  (3-3)  is  obtained  by  putting  180°  £  6_<360°,  180°  £  $  £  360° 

in  the  negative  of  one  of  the  equations  (3-36)  to  (3-39). 

Again  the  mathematical  evaluation  of  (3-36)  to  (3-39)  can  be  found 

» 

in  the  report  [5].  A  computer  program,  complete  with  operating  instruc¬ 
tions,  for  computing  the  transmission  through  a  rectangular  slot  in  a 
conducting  plane  is  given  in  Part  Two  of  [5], 

3-6  REPRESENTATIVE  COMPUTATIONS 


A  number  of  representative  computations  using  the  above  matrix  solu¬ 
tion  is  given  in  the  report  [5].  We  summarize  some  of  these  results  here. 


42 


The  first  computations  were  made  for  a  narrow  slot,  of  width  X/20 

in  the  y  direction  of  variable  length  L  in  the  x  direction.  The  far-zone 

quantity  plotted  was  the  transmission  cross  section,  given  by  (3-12), 

where  H  is  the  component  of  magnetic  field  being  considered.  We  use 
m 

the  notation: 

2  2 

T.  =  2vr  |h. I  in  the  y  =  0  plane 

0y  m1  0' 

(3-40) 

2  2 

r  =  27rr“)H  )  in  the  x  =  0  plane, 

xx  m 1  x 

We  choose  Ay  =  X/20  so  that  q  =  1  in  (3-27)  and  (3-28).  In  this  case, 

there  are  no  y  directed  magnetic  current  expansion  functions  (3-28).  Thus, 

the  matrix  solution  (2-5)  for  the  magnetic  current  M  is  an  x  directed 

vector.  As  a  result,  =  0  in  the  y  =  0  plane  and  =  0  in  the  x  =  0 

plane.  In  other  words,  for  the  case  being  considered,  the  components  of 

H  orthogonal  to  those  in  (3-40)  are  zero. 

Figure  3-3  shows  plots  of  xA  and  T  for  x-directed  slots  of  width 

0y  xx 

X/20  and  length  (a)  L  =  X/4,  (b)  L  =  X/2,  (c)  L  =  3X/4,  and  (d)  L  =  X.  In 
all  cases  the  excitation  was  due  to  a  plane  wave  normally  incident  on  the 
conducting  plane  with  the  magnetic  field  in  the  x  direction.  Note  the 
large  transmission  cross  section  for  L  «  X/2,  case  (b) ,  due  to  the  slot 
being  near  resonance.  The  plots  of  i  are  of  the  same  form  as  scattering 
cross  section  from  the  complementary  conducting  strips,  as  known  from 
Babinet's  principle. 

Figure  3-4  shows  plots  of  the  equivalent  magnetic  current  in  the 
aperture  region  for  the  same  slots.  Since  M  =  z.  x  E,  they  are  also  plots 
of  the  tangential  component  of  E  in  the  slots.  Again  note  the  large 
value  of  M  for  the  case  L  =  X/2,  which  is  near  resonance.  Note  also  that, 
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Fig.  3-4.  Magnitude  and  phase  of  j M/ E 1 j ,  where  M  is  the  x-directed 
magnetic  current  and  E*  is  the  incident  electric  field, 
for  the  same  slots  as  for  Fig.  3-3,  (a)  L  =  A/4,  (b)  L  =  A/2, 
(c)  L  =  3  / 4,  (d)  L  =  X.  Circles  denote  magnitude,  tri¬ 
angles  denote  phase. 
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for  short  slots  (L  3A/4),  the  M  is  almost  equiphasal  and  closely  approxi¬ 
mated  by  a  half  sine  wave. 

Next,  computations  were  made  to  test  the  rate  of  convergence  of 
the  solution  as  the  number  of  subsections  was  increased.  A  slot  of 
width  A/10  and  length  2A  was  chosen  for  the  study.  Again  the  excitation 
is  a  plane  wave  normally  incident  on  the  conducting  plane  with  the 
magnetic  field  in  the  x  direction.  Figure  3-5  shows  plots  of  r  and 
for  the  cases  (a)  39,  (b)  19,  (c)  9,  and  (d)  4  triangular  expan¬ 
sion  functions  respectively.  Note  that  the  patterns  (a)  and  (b)  are 
essentially  the  same,  and  pattern  (c)  is  only  slightly  different. 

They  differ  appreciably  from  (d),  which  results  from  only  4  expansion 
functions.  The  difference  in  the  solutions  as  the  number  of  expansion 
functions  is  decreased  is  better  illustrated  by  plots  of  M,  as  shown 
in  Fig.  3-6.  These  are  for  the  same  cases  as  for  Fig.  3-5.  It  can  be 
seen  clearly  how  the  computed  equivalent  current  in  the  slot  region 
changes  as  the  number  of  subsections  is  reduced.  As  a  rule  of  thumb, 
for  near-field  quantities  (such  as  M)  one  should  use  subareas  of 
length  A/10  or  less  and  for  far-field  quantities  (such  as  t)  length  1/5 
or  less. 


3-7  DISCUSSION 

The  computer  program,  given  in  Part  Two  of  [5],  is  written  ex¬ 
plicitly  for  rectangular  apertures,  but  the  formulas  are  valid  for  any 
aperture  composed  of  rectangular  subsections.  Other  apertures,  such  as 
L-shaped,  T-shaped,  square  O-shaped,  etc.,  could  be  treated  by  appro¬ 
priately  changing  the  computer  program.  Apertures  of  arbitrary  shape 


Transmission  cross  section  when  the  number  of  expansion  func¬ 
tions  is  (a)  19,  (b)  19,  (c)  9,  and  (d)  4.  Computations  are 
for  a  slot  of  length  21  in  the  x  direction  and  width  1/10  in 
the  y  direction.  Excitation  is  by  a  plane  wave  normally 
incident  on  the  conducting  plane  with  magnetic  field  in  the 
x  direction. 
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Fig.  3-6.  Magnitude  and  phase  of  |M/E  | ,  where  M  is  the  x-directed 
magnetic  current  and  is  the  incident  electric  field, 
when  the  number  of  expansion  functions  is  (a)  39,  (b)  19 
(c)  9,  and  (d)  4.  Circles  denote  magnitude,  triangles 
denote  phase.  Computations  are  for  the  same  slot  as 
for  Fig.  3-5. 
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could  be  treated  by  approximating  them  by  rectangular  subsections.  As 
with  all  moment  solutions,  the  size  of  the  apertures  which  can  be  treated 
depends  upon  the  size  of  the  matrix  which  can  be  computed  and  inverted. 

The  examples  indicate  that  the  rectangular  subsections  should  have  side 
lengths  not  greater  than  0.2  wavelengths  for  reasonable  accuracy. 

The  aperture  admittance  matrix  for  radiation  into  half-space  has 
application  to  any  problem  in  which  one  region  is  bounded  by  a  plane 
conductor,  as  shown  in  Chapter  II.  Hence,  it  can  be  used  for  a  waveguide- 
fed  aperture  in  a  ground  plane  ar.d  for  a  cavity-backed  aperture  in  a 
ground  plane.  These  problems  are  considered  in  the  next  two  chapters. 
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Chapter  IV 

WAVEGU IDE-FED  APERTURES 

4-1  GENERAL  THEORY 

Consider  now  a  uniform  waveguide  feeding  an  aperture  in  a  con¬ 
ducting  plane,  as  shown  in  Fig.  4-1.  In  general,  the  aperture  may  be 
of  different  size  and  shape  than  the  waveguide  cross  section.  The  halt- 
space  region  z >  0  is  the  same  as  in  the  previous  problem,  Fig.  3-1,  and 
the  analysis  of  the  preceding  chapter  applies.  An  analysis  of  the 
waveguide  region  is  given  below. 


Fig.  4-1.  Waveguide-fed  aperture  in  a  conducting  plane. 


Lot  the  excitation  of  the  waveguide  be  a  source  which  produces 
a  single  mode,  of  unit  amplitude,  incident  on  the  aperture.  This  mode 
(usually  the  dominant  mode)  is  denoted  by  the  index  o.  The  field  trans¬ 
verse  to  the  z-direction  can  then  be  expressed  in  modal  form  as 
(Sect.  8-1  of  [1]). 

-■)  z  y  z 

E=e  e  +  )  . .e  e. 

— t  — o  y  i  — i 


-Y  z  Y,z 

H  =  Y  e  u  v  e  -Vr.Y.e  u  *  e, 

— t  o  — z  o  “li  — z  — i 


(4-1) 


It  is  assumed  that  all  modes,  TE  and  TM,  are  included  in  the  summation. 
The  are  modal  propagation  constants 


Y  j  = 


j3.  =  Jk  A  -  (M\iV 


(4-2) 


(.  -i.  =  k.  A  -  (A  ./A) 2  A  >  At 

where  A^  is  the  ith  mode  cut-off  wavelength,  and  =  2ir/A^,  is  the  ith 
mode  cut-off  wavenumber.  The  Y.  are  the  modal  characteristic  admittances 


Y . 
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Yj/jwu 

jwc/Yi 


TE  modes 


TM  modes 


(4-3) 


is  the  reflection  coefficient  for  the  oth  mode,  and  1'  is  the. complex 
o  i 

amplitude  of  the  -z  traveling  component  of  the  ith  mode.  The  e^  are 
normalized  modal  vectors,  so  that  the  modal  orthogonality  relationships 


are 


jj  %  *  V"  =  < 

guide 


0 


i  *  j 
i  =  j 


(4-4) 


where  the  integration  is  over  the  waveguide  cross  section. 
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To  evaluate  the  aperture  admittance  (2-9)  for  the  waveguide  region, 
we  consider  a  single  expansion  function  on  the  z=0  plane  in  the  wave¬ 
guide  region.  The  tangential  field  produced  by  will  be  of  the  form 
(4-1),  except  that  there  is  no  incident  wave.  Hence,  this  field  is 
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where  the  superscript  wg  denotes  waveguide.  The  of  (4-9)  is  given 
by  the  second  equation  of  (4-5)  evaluated  at  z  =  0,  So  that 


YWg  =  l  A  .Y 
mn  u 


ni  l 


W 
— m 


apert . 


u  x  e.ds 
— z  —i 


(4-10) 


Now  define  the  constants 
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W 
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(4-11) 


which  are  similar  in  form  to  the  A  .  of  (4-8) .  The  elements  (4-10) 

m 

then  are  given  by 


YWg  =  J  A  .B  Y. 
mn  .  ni  mi  l 


(4-12) 


l 

Hence,  all  elements  of  the  waveguide  aperture  admittance  matrix  [Ywg]  are 

linear  combinations  of  the  modal  characteristic  admittance  Y  .  For 

Galerkin's  method,  W  =  M  and  A  ,  and  B  .  are  equal. 

— n  — n  ~ni  — ni 

We  next  evaluate  the  equivalent  magnetic  current  M,  given  by  (2-5). 
The  incident  field  is  given  by  the  first  term  on  the  right-hand  side  of 
(4-1).  When  the  aperture  is  covered  by  a  conductor,  the  waveguide  is 
terminated  by  a  conducting  plane.  According  to  image  theory,  the  tan¬ 
gential  magnetic  field  at  z  =  0  is  then  just  twice  the  incident  wave  or 


H1  =  2Y  u  x  e 
— t  o— z  — o 

This  is  the  used  in  (2-11)  to  evaluate  the  excitation  vector 
Hence,  the  components  of  the  excitation  vector  are 


(4-13) 
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The  total  aperture  admittance  matrix  is 

[Ya  +  Yb]  =  [YWg  +  Yhs]  (4-15') 

we  hs 

where  [Y  6]  is  the  waveguide  aperture  admittance  and  [Y  ]  is  the  half- 

-f 

space  aperture  admittance.  The  coefficient  matrix  V  is  given  by  (2-14) 
with  the  admittance  matrix  given  by  (4-15),  or 
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(4-20) 


(4-21) 


(4-22) 


The  parameter  of  most  Interest  Is  F  ,  the  reflection  coefficient  of 
the  incident  mode.  This  is  often  expressed  in  terms  of  an  admittance 


i  -  r 

=  - —  y 

tp  i  +  r  o 
o 


(4-23) 


which  is  the  equivalent  aperture  admittance  seen  by  the  incident 
mode. 

The  region  z  >  0  for  the  waveguide-fed  aperture  is  the  same  half¬ 
space  region  as  existed  in  the  previous  problem  of  an  aperture  in  a 
conducting  plane.  Hence,  evaluation  of  the  fields  in  terms  of  M  in  this 
region  is  done  in  the  same  way  as  in  Chapter  Ill.  For  example,  the  u^ 
component  of  the  far-zone  magnetic  field  at  a  point  r  is  given  by 


H  =  c  ^  ra  pm[YWg  +  Yhs]_1  I1  (4-24) 

m  4irr 

m 

which  is  (3-11)  with  the  term  [2Y^S]  replaced  by  (4-16).  The  excita¬ 
tion  vector  T *  has  elements  given  by  (4-14),  and  the  far-field  measure¬ 
ment  vector  P™  has  elements  given  by  (3-10).  The  power  gain  pattern  G  is 
the  ratio  of  the  radiation  intensity  in  a  given  direction  to  the  radi¬ 
ation  intensity  which  would  exist  if  the  total  power  Re(Pt)  were  radiated 
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uniformly  over  half  space,  or 


2Ttr2n|H  l2 

m  '  m ' 

Re(Pt) 


(4-25) 


Substituting  for  H  from  (4-24),  we  have 
m 


2  2 
oi  e  n 

8rrRe(Pt) 


|pm[YWg  +  Yhs]_1  I1 | 2 


(4-26) 


where  P  is  given  by  (2-27).  Note  that  this  gain  is  a  function  of 
the  II  component  measured,  as  well  as  direction  to  the  field  point. 


4-2  APPLICATION  TO  A  RECTANGULAR  WAVEGUIDE 

We  now  apply  the  general  theory  to  a  rectangular  waveguide  feeding 

a  rectangular  aperture  in  a  conducting  screen.  Figure  4-2  shows  the 

problem  to  be  considered  and  defines  the  coordinates  and  parameters  to 

be  used.  The  perfectly  conducting  plate  covers  the  entire  z  =  0  plane 

except  for  the  aperture  which  is  rectangular  in  shape  with  side  lengths 

L  Ax  and  L  Ay  in  the  x  and  y  directions  respectively.  L  and  L  are 
x  y  x  y 

positive  integers  and  ^  2.  The  aperture  is  fed  by  a  rectangular 
waveguide.  The  excitation  of  the  waveguide  is  a  source  which  produces 
one  mode,  of  unit  amplitude,  which  travels  toward  the  aperture. 

The  general  method  of  solution  discussed  in  Chapter  II  is  to  cover 
the  aperture  with  a  perfect  electric  conductor,  to  place  magnetic  cur¬ 
rent  sheets  +M  and  -M  respectively  on  the  left-hand  and  right-hand  sides 
of  this  conductor,  to  obtain  an  integral  equation  for  M  by  equating  the 
tangential  magnetic  fields  on  both  sides  of  this  conductor,  and  to  solve 
this  integral  equation  using  the  method  of  moments.  The  testing  func- 

i 

tions  are  the  same  as  the  expansion  functions  for  M  and  are  denoted  by  M.. 
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Fig.  4-2.  A  rectangular  waveguide  radiating  through  a  rectangular 

aperture  into  half-space  bounded  by  an  electric  conductor. 


Each  is  a  triangle  in  the  direction  of  current  flow  and  a  pulse  in 
the  direction  perpendicular  to  current  flow. 

Expression  (4-8)  for  A  requires  a  knowledge  of  the  expansion 
functions  and  waveguide  modes  e  .  The  set  of  expansion  functions 
is  split  into  a  set  of  x  directed  magnetic  currents  and  a  set  of 
y  directed  magnetic  currents  defined  by 
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where  T*(x)  and  T^(y)  are  triangle  functions  defined  by  (3-29)  and 
(3-30),  and  P^(x)  and  P^(y)  are  the  pulse  functions  defined  by 
(3-31)  and  (3-32).  An  evaluation  of  the  matrix  [YW®]  is  given  in 
detail  in  [2].  An  evaluation  of  the  matrix  [Y  ]  is  the  same  as 
that  for  the  rectangular  aperture,  given  in  Sect.  3-4  of  the  pre¬ 
ceding  chapter.  The  listing  and  documentation  of  a  computer  program 
for  the  problem  of  Fig.  4-2  is  given  in  Part  II  of  the  report  [2]. 


4-3  SAMPLE  COMPUTATIONS 

In  this  section  we  give  some  representative  computations  for 
the  aperture  of  Fig.  4-2.  Figure  4-3  shows  computed  results  for  a 
rectangular  waveguide  of  dimensions  A  by  A/2  radiating  into  half 
space  through  a  narrow  centered  rectangular  slot  of  dimensions  A  bv 
A/10,  that  is,  a  =  A  and  b  =  A/2.  Figure  4-3(a)  shows  the  x-component 
of  equivalent  magnetic  current,  which  is  also  equal  to  the  y-component 
of  tangential  E  field  in  the  slot.  No  y-component  of  magnetic  current 
was  obtained  because  only  one  pulse  in  y  was  used.  M  is  norma!  ' zed 
with  respect  to 

| e  |  dx  dy  (4-33) 

guide 

where  the  integral  is  over  the  waveguide  cross  section.  In  other  words, 
the  normalization  factor  is  the  root-mean-square  value  of  the  K  field 
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of  the  incident  wave.  The  phase  of  M  is  with  respect  to  that  of  this  E 
field  at  the  aperture.  All  computations  are  for  dominant  TE^  mode 
excitation.  Figure  4-3 (b)  shows  the  radiation  gain  patterns  in  the 
two  planes  x  =  0  and  y  =  0.  The  notation  Gg  denotes  the  gain  pat¬ 
tern  due  to  H.  in  the  y  =  0  plane.  The  notation  G  denotes  the 
0  xx 

gain  pattern  due  to  in  the  x  =  0  plane.  The  horizontal  axis  in 
Fig.  4-3 (b)  is  the  z  axis. 

Figure  4-4  shows  a  plot  of  the  equivalent  aperture  admittance 
Y  seen  by  the  dominant  mode  of  an  open-ended  square  waveguide  of 
width  a  radiating  into  half  space.  It  is  defined  by  (4-23),  where 


Fig.  4-4.  The  equivalent  aperture  admittance  seen  by  the  dominant  mode 
for  an  open-ended  square  waveguide  of  width  a  radiating  into 
half  space.  Our  computed  results  are  compared  to  those  cal¬ 
culated  and  measured  by  Cohen,  Crowley,  and  Levis  [3], 
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I’  is  the  reflection  coefficient  and  Y  is  the  characteristic  wave 
o  o 

impedance,  both  for  the  dominant  mode.  Our  computations  are  com¬ 
pared  to  some  previously  obtained  by  Cohen,  Crowley,  and  Levis  [3]. 
Also  shown  are  measured  values  reported  in  [3].  Additional  numerical 
computations  are  given  in  the  report  [2]. 
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Chapter  V 

CAVITY-RACKED  APERTURES 

5-1  THEORY  FOR  A  CYLINDRICAL  CAVITY 

An  important  special  case  of  the  cavity-backed  aperture  is  that 
where  the  cavity  is  a  finite  cylinder  of  arbitrary  cross  section.  The 
aperture  exists  at  one  end  of  the  cylinder,  and  the  other  end  is  com¬ 
pletely  covered  by  a  conductor.  The  cavity  can  then  be  viewed  as  a 
short-circuited  waveguide  and  waveguide  theory  applied.  Figure  5-1 
represents  a  typical  problem,  where  the  excitation  is  by  an  incident 
plane  wave  from  the  half -space  region. 


Fig.  5-1.  A  cavity-backed  aperture  where  the  cavity  is 
formed  by  a  shorted  waveguide. 
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The  half-space  part  of  the  problem  is  identical  to  that  treated 

in  Chapter  III.  For  a  unit  incident  plane  wave  of  polarization  u_^ 

and  propagation  vector  k  ,  the  excitation  vector  remains  that  of  (3-5). 

For  a  given  set  of  expansion  and  testing  functions,  the  aperture 

admittance  matrix  of  region  a  remains  the  same  as  that  for  Section  3-1, 
hs 

denoted  by  [Y  ] .  The  new  aspects  of  the  problem  are  those  of  deter¬ 
mining  the  aperture  admittance  matrix  for  the  cavity  region,  which 
we  denote  by  [Ycav] ,  and  the  measurement  vectors  I™  for  various  de¬ 
sired  field  quantities. 

The  field  in  the  cavity  region  can  be  expanded  in  terms  of  short- 
circuited  waveguide  modes.  For  a  cavity  depth  d,  the  transverse  to  z 
components  of  an  arbitrary  field  can  be  written  as 


E  =  -  I  « • 


sin  k . (d-z) 


— t  .  i  — i  sin  k.d 

l  l 


cos  k  (d-z) 

H  =  j  I  A.Y .  (u  x  e.)  - -  - 

— t  J  -  l  l  — z  —i  sin  k.d 

l  i 


(5-1) 


(5-2) 


Here  e.  are  the  normalized  modal  electric  field  vectors  (discussed  in 

—i 

Chapter  IV,  Y^  are  the  modal  characteristic  admittances,  and  k^  are  the 
modal  wave  numbers.  The  term  sin  k^d  is  placed  in  the  denominator  for 
later  convenience.  The  summations  in  (5-1)  and  (5-2)  are  assumed  to  be 
over  all  modes,  both  TE  and  TM  if  necessary.  In  particular,  the  modal 
wave  numbers  are 


5  --  ) 

ic 


A  >  A. 

ic 


(5-3) 
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where  AJ  =  2Ti/k,  is  the  ith  mode  out-off  wavelength.  The  modal 
ic  ic 

characteristic  admittances  are 


,  TE  modes 

^toe/k^  ,  TM  modes 


(5-4) 


The  intrinsic  wavenumber  k  may  be  complex  if  u  and/or  £  are  complex 
to  account  for  a  lossy  medium  in  the  cavity.  The  e^  are  normalized 
modal  vectors,  so  the  modal  orthogonality  relationships  are 

0  ,  i  *  j 


e .  •  e .  ds  = 

-i  -j 


(5-5) 


1  , 


Here  the  surface  of  integration  S  is  the  cavity  (waveguide)  cross  section. 

To  evaluate  the  aperture  admittance  (2-10)  in  the  cavity  region,  we 
consider  a  single  expansion  function  on  the  z=0  plane  inside  the 
cavity.  The  only  source  in  the  cavity  region  is  M^,  hence  the  tangential 
field  produced  by  is  of  the  form  of  (5-1)  and  (5-2).  Using  the  bound¬ 

ary  condition  M  =  -  u  x  E  and  specializing  (5-1)  to  z  =  0,  we  have 


M  =  -  u  x  K 
— n  — z  — t 


z  =  0 


=  I 


A  .u  x  e . 
ni—z  —i 


(5-6) 


Here  the  additional  subscript  n  is  placed  on  A  ,  to  denote  that  it  is 

ni 

that  due  to  M  .  Multiplying  each  side  of  (5-6)  by  u  x'  <2.  and  inte- 
~ n  z  j 

grating  over  the  cross  section  S,  we  obtain 


M  •  u  x  e.ds  =  j  A  , 

JJ  -n  -z  -j  L  ni  Jj 
S  S 

By  orthogonality  (5-5),  all  terms  of  the  summation  are  zero  except  the 


<u 


X  ^i} 


e  )ds 


(5-7) 


i“j  term.  Hence 
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A  =  M  •  u  x  e.ds  (5-8) 

nl  j  j  n  z  1 

apert . 


We  have  replaced  the  integral  over  S  by  one  over  the  aperture  because 
exists  only  in  the  aperture  region. 

The  elements  of  the  aperture  admittance  matrix  (2-10)  for  the 
cavity  region  are  now  given  by 


ycav 

mn 


W  •  HT (M  ) ds 


(5-9) 


apert. 


The  H  of  (5-9)  is  given  by  (5-2)  evaluated  at  z=0,  so  (5-9)  becomes 


YCaV  =  -  j  l  A  .Y.  cot(k.d) 
mn  J  j  ni  i  1 


W  •  u  x  e.ds 
— m  — z  — i 


(5-10) 


apert . 


We  now  define  the  constants 


B 


mi 


apert. 


x  e 
% 


ds 


(5-11) 


which  are  similar  in  form  to  A  ^  given  by  (5-8) .  Then  the  admittances 
(5-10)  are  gi\en  by 

Ycav  =  -  j  )"  A  B  Y  cot(k.d)  (5-12 

mn  “  ni  mi  i  i 

C3  V 

Hence,  the  elements  of  [Y  ]  are  linear  combinations  of  the  input  wave¬ 
guide  admittances  for  each  short-circuited  waveguide  mode. 

Some  specific  quantities  of  interest  in  the  solution  are  (a)  the 
equivalent  magnetic  current  M,  or  tangential  15  in  the  aperture,  (b)  the 
amplitude  of  some  specific  mode,  and  (c)  the  electric  field  intensity 
at  some  point  in  the  cavity.  Once  M  is  obtained,  then  quantities  (b) 
and  (c)  are  easily  obtained. 
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We  have  the  equivalent  magnetic  current  given  by  (2-5)  where  V 

n 
— y 

are  the  elements  of  V,  obtained  from  (2-14)  specialized  to  the  present 
problem.  This  result  is 


V  =  [Yhs  +  Y^V1  P1 


(5-13) 


where  the  elements  of  P*  are  given  by  (3-5).  The  modal  amplitudes 
are  obtained  from  (5-6)  as  it  applies  to  the  total  magnetic  current  in 
the  cavity,  which  is  -M.  The  result  is 


-M  =  ^  A .  u  x  e . 


(5-14) 


Again  we  multiply  each  side  by  u  x  ,  and  integrate  over  the  cross 

2  J 

section  S,  as  in  (5-7),  to  obtain 


Ai 


(-M)  *  u  x  eds 


apert. 

Substituting  for  M  from  (2-5),  and  using  (5-13),  we  have 


(5-15) 


At  =  -  AV  =  -AlYhs  +  Y^V1  P1 


(5-16) 


where  A  is  the  row  vector  with  elements  given  by  (5-8).  Hence,  the 

-+ 

measurement  vector  for  obtaining  mode  amplitudes  is  the  vector  A. 

Finally,  to  obtain  the  field  at  any  point  in  the  cavity,  we 
use  (5-1)  with  the  A^  given  by  (5-16).  At  a  frequency  near  resonance, 
that  is,  where  sin  k^d  ~  0  for  some  j,  a  good  approximation  to  the 
field  (5-1)  is  given  by  only  the  jth  term  of  (5-1),  or 


-  A ,  e . 


sin  kj (d-z) 


— t  j— j  sin  k.d 

J  J 


(5-i, 
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If  the  cavity  is  truly  loss  free,  then  (5-17)  seems  to  predict  an 

infinite  field  at  resonance.  This,  however,  is  not  correct,  since 

A,  +  0  as  sin  k,d  -*■  0  giving  a  finite  value  for  E  at  resonance, 
j  j  “t 

The  excitation  of  the  resonant  mode  must  then  be  determined  from  the 
orthogonality  condition 


M  •  u  x  e.  ds  =  0  (5-18) 

J  J  -  -z  -3 

apert . 

However,  the  theory  for  this  special  case  will  not  be  considered  here, 
since  it  requires  considerable  modification  of  the  formulas  of  this 
section. 


5-2.  APERTURE  BACKED  BY  AN  ARBITRARY  CAVITY 

When  the  cavity  backing  the  aperture  is  of  arbitrary  shape,  there 
are  several  ways  of  calculating  the  field  in  the  cavity.  One  way  is 
completely  modal.  This  is  the  method  used  in  this  chapter.  Another 
way  is  a  nonmodal  approach,  similar  to  that  used  in  [2],  That  method, 
however,  will  not  be  considered  here. 

Figure  5-2  represents  the  cavity  part  of  a  problem  of  the  type  of 


A 

n 


Fig.  5-2.  Equivalence  for  cavities  of  arbitrary  shape. 
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Fig.  5-1,  except  that  the  cavity  is  now  of  arbitrary  shape.  It  is 
desired  to  calculate  the  field  _E  and  H  in  the  cavity  region  R,  com¬ 
pletely  enclosed  by  a  conducting  surface  S,  and  excited  by  the 
equivalent  current  M  over  the  aperture  surface.  In  general,  the 

magnetic  field  in  R  can  be  represented  as  the  sum  of  a  curl-free 

c  f 

(irrotational)  part  H  plus  a  divergence-free  (solenoidal)  part 
Hd^,  that  is 

H  =  Hcf  +  Hdf  (5-19) 


The  division  (5-19)  is  not  necessarily  unique,  since  part  of  the  field 
may  be  both  curl-free  and  divergence-free.  For  simplicity,  we  re¬ 
strict  consideration  to  simply-bounded  and  simply-connected  cavities. 

The  modal  representation  of  the  magnetic  field  in  a  cavity 
bounded  by  a  perfect  electric  conductor  is  available  in  textbooks 
[3],  [4].  The  use  of  this  representation  for  aperture  problems  is 
summarized  well  in  [5].  The  divergence-free  part  of  the  field  can  be 
represented  in  terms  of  resonant  cavity  modes  as 


af  IL(£)  ff 

H  =  ju>e  l  ^—2  (- 

lk"ki  aper t . 


(-M)  •  ds 


(5-20) 


The  are  solutions  to  the  eigenvalue  problem 


V  x  V  x  Hi  =  k^H£  in  R 


(5-21) 


nxyxHi=0  onS  (5-22) 

2  2 

The  eigenvalues  ^  are  real,  and  only  k1  >  0  are  used  in  (5-20).  The 
eigenvectors  IT  are  also  real,  and  orthonormal  according  to 


r 

\ 


\ 
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'ff 

J. 

R 


Mi 


H.  ds  = 


1 


i  *  i 
i  =  j 


(5-23) 


The  H  defined  above  are  the  usual  resonant  modes  of  time-harmonic 
fields  in  the  cavity. 

The  curl-free  part  of  the  field  can  be  expressed  in  terms  of 
"static  modes"  as 


,cf 


~  y  G.(r)  [[  (-M)  .  G .  ds 
jojy  £  -i  -  j  “  -i 


apert . 


(5-24) 


where 


Si- 


(5-25) 


-v\  = 


in  R 


(5-26) 


8n 


=  0 


on  S 


(5-27) 


The  modal  vectors  satisfy  (5-21)  and  (5-22)  for  *  0,  that  is, 

2 

for  zero  eigenvalue.  Hence,  we  can  think  of  ^  -  0  as  being  an  eigen¬ 
value  of  infinite  degeneracy.  The  G.  are  real  and  orthonormal  accord- 

— i 


ing  to 


r 


G  •  G  ds 
“i  -j 


i  *  J 


i  =  i 


(5-28) 


Also,  since  is  orthogonal  to  HC^  in  general,  we  have 


G,  •  H,  ds  =  0 


(5-29) 


for  all  i  and  j.  The  total  field  in  R  is  simply  the  sum  of  (5-20) 
and  (5-24)  according  to  (5-19). 


69 


There  is  an  alternative  form  for  (5-24)  which  is  sometimes  con¬ 
venient.  This  is 


c  f 

H  =  -  Vi|> 


(5-30) 


where 


'P  =  y  ^i(-)  |  ds 


(5-31) 

i. 

apert. 

Here  m  is  the  magnetic  charge  associated  with  M  according  to  the  equa¬ 
tion  of  continuity 


m  =  - —  V  •  M 

jo)  -  - 


(5-32) 


The  ( if/.  are  still  solutions  to  the  eigenvalue  problem  (5-26)  and  (5-27), 
and  orthonormal  according  to  (5-28)  where  G^  is  given  by  (5-25).  The 
derivation  of  (5-31)  from  (5-24)  involves  substitution  for  C.  from 
(5-25)  and  application  of  the  divergence  theorem  to  the  integral  over  R. 
As  shorthand  notation,  we  can  define  the  set  of  modes  {F.;  lo  be 

— l 

{Gl5  j^},  that  is,  to  include  all  modes.  Then  the  total  magnetic  field 
can  be  written  as 


F.(r) 

H  =  jwG  l 


2  2 
i  k  -k, 

i  apert 


where 


F  ♦  F.  ds  = 
-i  -1 


( -M)  •  F  d  s 


i  *  j 


(5-33) 


(5-34) 


This  is  the  form  used  in  Sec.  8-13  of  []].  However,  we  must  romember 
that  there  are  an  infinite  set  of  modes  associated  with  the  eigenvalue 

kj  -  0. 

To  obtain  the  aperture  admittance  elements  for  the  cavity  region 


we  U3e  (5-9)  where  H  is  now  given  by  (5-33)  with  M  replaced  bv  M  ,  or 
—  n 
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the  sum  of  (5-20)  and  (5-24)  with  M  replaced  by  M^.  The  result  is 
,cav  1 


b\  b  , 
mi  ni 


Y  =  - -  J  a',  a  +  jco£  T  ,  „ 

mn  mi  ni  i  k2  -  k2 


(5-35) 


where 


a  . 
ni 

= 

J 

‘  M 
— n 

a*  . 

1  ni 

apert . 

W 

L-n 

G .  ds 
— i 


(5-36) 


b  . 
ni 

. 

• 

j 

M 
— n 

b', 
ni  . 

apert . 

W 

L~*n 

•  H,  ds 


(5-37) 


Once  again  we  have  the  possibility  that,  for  a  loss-free  cavity,  one 
or  more  terms  of  (5-35)  can  become  infinite  at  resonance,  that  is, 
when  k  =  k. .  In  the  case  of  lossy  cavities  this  does  not  occur,  although 
numerical  problems  may  arise  if  the  cavity  is  only  slightly  lossy. 

The  magnetic  field  produced  by  the  nth  expansion  function  is  given 
by  substituting' into  the  sum  of  (5-20)  and  (5-24),  or 


h(m  )  =  ~  y  a  ,g,  +  jwe  y 

—  ~n  igju  i  ni— i  L, 


b  4  H. 
ni  —l 

2  2 

i  k  -  k^ 


(5-38) 


where  a  .  and  b  .  are  given  by  (5-36)  and  (5-37).  The  electric  field 
ni  ni 

is  obtained  from  (5-38)  according  to  jweE  =  V  x  H,  resulting  in 


E(M  )  =  l 


b  .  V  *  H. 
ni  —  — i 


2  2 
i  k  -  kt 

Note  that  the  first  summation  of  (5-38)  vanishes  when  (5-39)  is 


(5-39) 


derived,  since  7  *  =  0.  Both  E  and  H  are  linearly  related  to 
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magnetic  current.  Remembering  that  the  source  in  region  b  of  the  original 
problem  is  -M,  we  can  use  the  superposition  (2-5)  to  obtain  the  total  mag¬ 
netic  and  electric  fields  in  the  cavity  as 


l  V  H(M  ) 

L  n - n 

n 


(5-40) 


E  =  -  T  V  E(M  ) 

—  u  n - n 

n 


(5-41) 


Here  the  elements  of  V  are  obtained  from  (2-14)  in  general.  For  the  par- 
ticular  case  of  a  cavity-backed  aperture  in  a  conducting  plane,  V  is 
given  by  (5-13)  where  [Ycav]  now  has  the  elements  (5-35).  In  (5-40),  the 
coefficient  of  the  jth  static  mode  Gj  is 

a.  =  "t-—  Y  V  a  .  =  a.V  (5-42) 

3  3<jy  “  n  nj  jugi  j 

where  a,  is  the  row  vector  with  elements  a  ,  given  by  (5-36).  Using 
J  nj 

(5-13)  for  V,  we  have  for  the  amplitude  of  the  jth  static  mode 


a  *  a .  [Yhs  +  yc3v]-1  ^ 

j  i^y  i 


(5-43) 


Similarly,  the  coefficient  of  the  jth  resonator  mode  in  (5-40)  is 


b.  -4S-IV  b  ,  .  JKEyVf 
J  k2-k2  „  “  "3  w2-k2  3 


(5-44) 


k  -k: 


where  b.  is  the  row  vector  with  elements  b  ,  given  by  (5-37).  Again 
J  nj 

using  (5-13)  for  V,  we  have  for  the  amplitude  of  the  jth  resonator  mode 


b,  ■  b.iv1" .  Y«vr1  ;* 

1  k  j 
»/,2  .  2. 


(5-45) 


Hence,  a^/juju  and  -ju)ebj/(k  -k^)  are  the  measurement  vectors  for  deter¬ 
mining  a  mode  amplitude. 


72 


5-3  DISCUSSION 

A  general  formulation  for  a  cavity-backed  aperture  in  an  infinite 
conducting  plane  has  been  given  in  terms  of  generalized  network  parameters. 
The  approach  is  also  valid  for  a  cavity-backed  aperture  in  a  conducting 
body  of  arbitrary  shape,  but  then  the  calculation  of  the  excitation  and/or 
the  measurement  vectors  is  more  difficult. 

Two  interesting  specializations  of  the  general  theory  are  (a)  the 
low  frequency  case  where  cavity  dimensions  are  small  compared  to  wavelength 
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and  (b)  the  resonant  case,  particularly  in  the  vicinity  of  the  first 
resonance.  It  is  hoped  that  the  general  theory  can  be  simplified  for 
these  cases,  allowing  us  to  derive  relatively  simple  equivalent  circuits. 
These  topics  require  further  investigation. 
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Chapter  VI 

ELECTROMAGNETIC  TRANSMISSION  THROUGH  NARROW  SLOTS 
IN  THICK  CONDUCTING  SCREENS 


6-1  INTRODUCTION 

The  first  accurate  treatment  of  electromagnetic  coupling  through 
small  holes  was  for  zero-thickness  conductors,  and  used  the  concept  of 
aperture  polarizability  [1],  [2].  This  theory  has  found  extensive  appli¬ 
cation  in  the  literature  (see  [3]  for  a  bibliography).  When  the  aperture 
becomes  larger  and/or  the  conductor  has  thickness,  the  polarizability 
concept  becomes  inadequate  and  the  generalized  admittance  concept  should 
be  used.  In  this  chapter  we  apply  the  admittance  concept  to  narrow, 
infinitely  long  slots  (sometimes  called  slits)  in  a  conducting  plane  of 
finite  thickness.  The  concepts  we  use  are  general,  applying  also  to 
three-dimensional  problems,  but  we  here  consider  only  the  two-dimensional 
case . 

The  problem  of  a  slot  in  a  thick  conducting  screen  has  been  con¬ 
sidered  by  several  methods  (see  [4]  for  a  bibliography).  It  was  treated 
by  the  generalized  admittance  concept  for  the  TE  (transverse  electric 
to  the  slot  axis)  case  in  [4].  A  similar  solution  for  the  TM  (transverse 
magnetic  to  the  slot  axis)  case  can  be  found  in  [5].  For  narrow  slots 
in  thick  conductors,  only  the  transmission  line  mode,  which  is  n  TE  mode, 
can  propagate  through  the  slot  region.  All  higher  order  TE  modes,  and 
all  TM  modes,  are  cut  off.  Hence,  we  restrict  explicit  consideration 
to  the  TF.  case.  The  solution  used  is  basically  a  simplification  of  the 
general  solution  [4]  for  the  narrow  slot.  This  results  in  relatively 
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simple  formulas,  and  a  simple  equivalent  circuit,  for  narrow  slots  in 
conducting  screens  of  finite  thickness. 


6-2  FORMULATION  OF  THE  PROBLEM 


The  problem  to  be  considered  is  shown  in  Fig.  6-1,  which  represents 
the  cross  section  of  a  conducting  screen  of  thickness  d  in  which  a  slot 
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Fig.  6-1.  A  uniform  slot  of  width  w  in  a  perfectly  conducting 
screen  of  thickness  d. 
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of  uniform  width  w  is  cut.  The  left-hand  half  space  (x  <  0)  is  called 
region  a,  the  uniform  slot  region  (0  <  x  <  d,  0  <  y  <  w)  is  called 
region  b,  and  the  right-hand  half  space  (x  >  d)  is  called  region  c.  The 
boundary  common  to  regions  a  and  b  is  called  the  aperture  A^.  The 
boundary  common  to  regions  b  and  c  is  called  the  aperture  .  Regions 
a,  b,  and  c  are  each  filled  with  homogeneous  media  of  constitutive 
parameters  (y  ,  £  ),  (y  ,  £,),  and  (y  ,  e  )  respectively.  Each  y  and  £ 

3  3  D  D  CC 

can  be  considered  complex  to  account  for  dissipation.  The  excitation 
is  due  to  known  sources  J_^  and  M'*'  in  region  a.  It  is  desired  to  obtain 
a  solution  for  the  field  in  each  region,  and  for  the  power  transmitted 
into  region  c. 

As  described  in  Chapter  II,  the  equivalence  principle  can  be  used  to 
divide  the  original  problem  into  three  equivalent  parts,  as  shown  in 
Fig.  6-2.  In  Fig.  6-2a,  we  have  the  original  sources  J1,  M1,  plus  the 
equivalent  magnetic  current  M^,  where 

M  =  x  x  E  (6-1) 

over  the  aperture  region  A^,  all  radiating  in  the  presence  of  a  complete 
conductor  (aperture  shorted).  In  (6-1),  x  is  the  x-directed  unit  vector, 
normal  to  A-^,  and  E  is  the  electric  field  in  the  aperture  A^  in  the  original 
problem.  In  Fig.  6-2b,  we  have  the  equivalent  magnetic  currents  -bL  ,  given 
by  (6-1),  over  A^,  and  -M^,  given  by 

M2  =  -  x  x  E  (6-2) 

over  the  aperture  region  A^,  all  radiating  in  the  presence  of  a  conductor 
completely  enclosing  the  rectangular  region  b  (both  apertures  shorted). 
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(c)  Equivalence  for  region  c. 


6-2.  The  problem  divided  into  three  equivalent  parts. 
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In  (6-2),  the  unit  vector  x  is  normal  to  A^,  and  II  is  the  electric  field 
in  the  aperture  A9  in  the  original  problem.  Finally,  in  Fig.  6-2c,  we 
have  the  equivalent  magnetic  current  given  by  (6-2),  over  the  aper¬ 
ture  region  A9 ,  radiating  in  the  presence  of  a  complete  conducting  plane 
(aperture  A^  shorted). 

The  use  of  in  Fig.  6-2a  and  -M^  in  Fig.  6-2b  ensures  the  conti¬ 
nuity  of  the  tangential  components  of  across  the  aperture  A^.  The  use  of 
-Mj  in  Fig.  6-2b  and  in  Fig.  6-2c  ensures  the  continuity  of  the  tan¬ 
gential  components  of  II  across  the  aperture  A^.  The  remaining  boundary 
conditions  to  be  enforced  are  the  continuity  of  tangential  components  of 
H  across  A^  and  A^ .  By  the  method  used  in  [A],  we  obtain  the  following 
two  operator  equations: 


-t(-l) 

-  -  h“ 

over  A^ 

(6-3) 

-t(-l) 

-  i£<«2> 

-  i£(H2>  -  o 

over  A^ 

(6—4) 

Here  h£(M^)  is  the  operator  which  gives  the  tangential  component  of  11  due 

SC 

to  the  current  M  radiating  in  region  p,  with  all  apertures  shorted.  H 
-q  — t 

is  the  tangential  component  of  11  due  to  the  impressed  sources  J  ,  M  in 

SC 

region  a  with  aperture  A^  shorted.  In  terms  of  the  current  J_  which 


exists  on  the  shorted  aperture  A^  in  Fig.  6-2a, 


H*C  =  x  x  jsc  (6-5) 

Sc 

In  Chapter  II,  H  was  called  the  field  due  to  the  impressed  sources 
and  denoted  by  H1.  Others  [6]  have  called  it  the  generator  field  and 
denoted  it  by  H  . 
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A  convenient  way  to  reduce  the  operator  equations  (6-3)  and  (6-4) 
to  matrix  equations  is  the  method  o£  moments  [7].  The  procedure  can 
be  summarized  as  follows.  Define  sets  of  expansion  functions  {M  }  in 
and  {l^  }  in  A ^ >  and  express  the  equivalent  magnetic  currents  as 


M  =  l  V,  M. 
—1  '-.In  —In 

n=l 


N2 

—2  '  J,  V2„  — 2n 
n=l 


(6-6) 


(6-7) 


where  V,  and  V.  are  coefficients  to  be  determined. 
In  2n 

product  for  each  slot  as 


<A,  B>q 


_B  dy 


Define  a  symmetric 


(6-8) 


where  q  ■  1  or  2.  The  integrand  in  (6-8)  is  evaluated  at  x  *  0  when 

q  =  1  and  at  x  =  d  when  q  =  2.  Define  sets  of  testing  functions 

(w,  }  in  A,  and  {w  )  in  A  .  Substitute  (6-6)  and  (6-7)  into  (6-3) 

— lm  1  —2m  2 

and  (6-4),  and  test  the  resultant  equations  with  each  W^,  m=l,2, . . .  ,N^, 

and  W.  ,  m=l,2, . . . ,N„.  The  result  is 

— zm  z 

IYUl7l  +  <YUi7l  +  [YUi52  ■  <6-9) 


where 


|Y21)51  +  +  IY22l52  ’  3 


[yP  ]  =  {-  <w  ,  hJ(m  )>  ]  xN 

qr  qm  t  rn  q  x 


Vt  ■  'V»  X  1 

r 


11  ■  V  x  i 


(6-10) 

(6-11) 

(6-12) 


(6-13) 


The  matrices  [YP  ]  are  called  the  generalized  admittances,  the  vectors 

-v  -►i 

are  called  the  generalized  voltages,  and  the  vector  I  is  called  the 

generalized  source  current.  A  solution  to  the  problem  is  obtained  by 
solving  the  matrix  equations  (6-9)  and  (6-10)  for  and  V^,  which  deter¬ 
mine  the  magnetic  currents  by  (6-6)  and  (6-7).  Once  the  equivalent 
magnetic  currents  are  known,  the  fields  in  each  region  can  be  obtained 
from  the  equivalent  problems  of  Fig.  6-2. 


6-3  SPECIALIZATION  TO  NARROW  SLOTS 

As  noted  in  the  introduction,  rigorous  solutions  to  the  general 
problem  can  be  found  in  the  literature.  We  here  consider  only  an 
approximation  which,  as  we  shall  see,  gives  highly  accurate  solutions 
for  narrow  slots.  We  consider  only  the  case  of  TE  excitation  because, 
as  noted  earlier,  all  TM  modes  in  a  narrow  slot  region  are  cut  off. 

The  approximate  solution  is  basically  a  one-term  moment  solution 
to  the  general  problem.  For  a  testing  function  in  each  aperture,  we 
choose  the  constant  (integrated  value  unity) 


-11  ‘  -21  '  -  ’  -  » 


(6-14) 


This  is  the  H  field  variation  of  the  transmission  line  mode  in  region  b, 
and  is  orthogonal  to  all  higher-order  waveguide  modes.  For  now,  the 
expansion  function  in  each  aperture 

-11  =  M21  =  M  =  i  f(y)  (6-15) 

will  be  left  undefined  as  to  its  functional  form  f(y).  However,  so 
that  M  excites  the  transmission  line  mode  the  same  regardless  of  its 


functional  form,  we  require 
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f (y)  dy  =  1 


(6-16) 


for  all  "trial  functions"  f(y).  In  other  words,  the  net  magnetic 
current  is  unity. 

The  general  network  equations  (6-9)  and  (6-10)  now  reduce  to  scalar 


equations 


YaV  4-  Yb  V  +  Yb  V  =  I1 
1  11  1  12  2 


(6-17) 


Yb  V  +  Yb  V  +  YCV  =  0 
21  1  22  2  2 


(6-18) 


where 


Y3  =  —  z  •  H*(M)  dy 
w  —  — t  — 


(6-19) 


is  the  aperture  admittance  of  looking  into  region  a,  and 


YC  =  — 


£  •  H  (M)  dy 


(6-20) 


is  the  aperture  admittance  of  looking  into  region  c.  The  matrix  of 
admittances 


Yb  Yb 

11  12 


Yb  Yb 

21  22 


-jY^cot  k^d  -jY^csc  k^d 


-jYQCsc  k^d  -jY^cot  k^d 


(6-21) 


is  the  two-port  admittance  matrix  for  a  parallel-plate  transmission  line 


of  length  d,  characteristic  admittance  Y  =  1/wn,  ,  and  propagation  constant 

o  b 

jk^,  where  and  are  the  intrinsic  impedance  and  wave  number  of  region  h. 
The  impressed  current  is 
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(6-22) 


which  is  the  average  surface  density  of  electric  current  over  the 
short-circuited  aperture.  The  approximations  to  the  equivalent 
magnetic  currents  in  the  apertures  and  are  now 

=  V  M  and  M2  =■  V2M  (6-23) 


respectively.  The  equivalent  circuit  for  this  approximation  to  the  narrow 
slot  is  shown  in  Fig.  6-3. 

It  is  well-known  that  the  aperture  admittance  of  a  thin  slot  open¬ 
ing  into  half  space  is  insensitive  to  small  variations  in  the  tangential 
electric  field  in  the  slot  about  its  true  value  [8],  We  here  consider 
three  cases,  all  of  which  result  in  an  aperture  admittance  of  the  form 

Ya  =  -y-  [it  -  2j  ln(C  k  w)]  (6-24) 

a  a 

for  region  a,  and  of  the  same  form  with  all  a's  replaced  by  c's  for 


Fig.  6-3.  Equivalent  circuit  for  a  narrow  slot  in  a  thick 
conducting  screen. 
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region  c.  In  (6-24),  n  ,  X  ,  and  k  are  the  intrinsic  impedance,  wave- 

3  3  3 

length,  and  wave  number  of  medium  a,  respectively,  and  C  is  a  constant 
depending  on  the  choice  of  f(y).  If  the  tangential  electric  field  in  the 
slot  is  assumed  constant,  then  [8] 

C  =  y/2e  '  ~  0.1987  (6-25) 

Here  y  =  1.781...  and  e  =  2.718...  .  If  the  tangential  electric  field 
is  assumed  to  be  the  quasi-static  solution  for  a  right-angle  flange, 
then  [9] 

C  =  Y/rre  ~  0.2086  (6-26) 

Finally,  if  the  tangential  electric  field  is  assumed  to  be  the  quasi¬ 
static  solution  for  a  slot  in  a  zero-thickness  screen,  then  [10] 

C  -  y/8  ^  0.2226  (6-27) 


Since  C  is  in  a  logarithmic  term  in  (6-24),  it  makes  little  difference 
which  value  is  chosen.  However,  so  that  our  solution  remains  strictly 
valid  as  the  thickness  d  -*•  0,  we  choose  (6-27). 

A  parameter  of  interest  is  the  transfer  admittance 


Y 


12 


(6-28) 


which  allows  one  to  calculate  the  strength  of  =  V^M  given  the  excita¬ 
tion  I  of  (6-22).  From  transmission  line  theory  and  circuit  theory  applied 
to  the  equivalent  circuit  of  Fig.  6-3,  we  obtain 

3  C 

Y12  *  -(Ya  +  YC)cos  l^d  -  J(Yq  +  Y-~)sin  k^d  (6-29) 

O  ) 

I 

! 
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The  power  transmitted  through  the  aperture  is  equal  to  that  dissipated 

c 

in  Y  of  the  equivalent  circuit,  that  is 


P 


trans 


V2|2  Re(YC) 


(6-30) 


In  terms  of  the  transfer  admittance,  this  is 


P 

trans 


Re(Yc) 


(6-31) 


where  I*  is  given  by  (6-22),  by  (6-29),  and  YC  by  (6-24)  with  a's 
replaced  by  c's. 

Now  consider  the  case  of  a  normally  incident  plane  wave 

.  -jk  x 

Hi0  =  H  e  3 
z  o 


(6-32) 


Here  the  superscript  io  denotes  a  wave  In  media  a  of  infinite  extent  (no 
conductor).  The  power  incident  on  the  aperture  (per  unit  length  in  the 
z  direction)  is  then 


P 


inc 


nJH 


w 


(6-33) 


The  short-circuit  magnetic  field  is  twice  the  incident  field  on  the  con¬ 
ductor,  or 

HSC  =  2H  (6-34) 

z  o 

Hence,  the  impressed  current  is,  by  (6-22), 

I1  =  2H  (6-35) 

o 

The  transmission  coefficient  of  the  slot  is  defined  as 


T 


P 


P 


trans 

inc 


(6-36) 


J 
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Now,  substituting  from  (6-31),  (6-33),  and  (6-35)  into  (6-36),  we  have 


4  Re(Y  ) 


(6-37) 


If  the  plane  wave  is  incident  at  some  angle  0.  in  the  x-y  plane  measured 

inc 

from  the  x  axis,  then  the  left-hand  side  of  (6-37)  becomes  T  cos  0. 

me 

The  transmission  coefficient  (6-37)  depends  on  screen  thickness  d 
only  through  the  parameter  Y  2»  In  particular,  it  will  be  maximum  when 
| Y^2 I  Is  minimum.  Consider  the  case  of  region  b  lossless,  and  denote  the 
aperture  admittances  by 


„a  _a  .  ,_a  .  ,,c  _c  ,  .  „c 

Y  =  G  +  jB  and  Y  =  G  +  jB 


(6-38) 


We  can  now  write  the  real  and  imaginary  parts  of  Y^  as  given  by  (6-29) 


Re(Y12)  =-(Ga  +  GC)cos  k^d  +  (G3BC  +  GCBa)sin  kfad 

o 


(6-39) 


gV  -  bV, 


Im(Y12)  =  — (B  +  B  )cos  k^d  -  (Yq  +  ------)sin  kbd 


(6-40) 


For  narrow  slots,  we  see  from  (6-24)  that  Ba  »  G3  and  BC  >>  GC.  Hence, 
the  coefficients  of  the  trigonometric  terms  in  (6-40)  are  much  larger 
than  those  in  (6-39),  and  we  can  minimize  lYb2|  by  setting  Im(Yj2)  =  0. 
As  w  -*•  0,  we  retain  only  the  dominant  Y^  term  in  the  coefficient  of 
sin  k^d  in  (6-40).  Then  Im(Y^2)  =  0  when 


,  J  _  ,Ba  +  BC, 
tan  kbd  ~  -  ( - - - ) 

o 

which  is  the  condition  for  "slot  resonance."  Since  Ba  +  Bl  ■  Y  , 

o 

a  first  approximation  to  resonance  is  l^d  nu ,  or 


(6-41) 
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d 

res 


-  \ 

2  V 


n  =  1,2,3, .. . 


(6-42) 


Here  the  subscript  "res"  denotes  "at  resonance."  Actually,  the  right- 
hand  side  of  (6-41)  is  a  small  negative  number,  hence  resonance  occurs 
when  d  is  slightly  less  than  an  integer  number  of  half  wavelengths  in 
region  b.  Assuming  that  tan  k^d  varies  linearly  in  the  vicinity  of 
each  zero,  we  have 


d 

res 


(n 


BS  +  BC.  \ 
ttY  ;  2 


(6-43) 


o 

where  n  =  1,2,3,...  .  This  is  a  better  approximation  to  the  resonant 
thickness  than  is  (6-42). 

We  next  wish  to  obtain  the  transmission  coefficient  at  resonance. 
Now  cos  k,  d  ~  (-l)n  and  sin 


k^d  is  small,  hence  from  (6-39)  we  have 


At  resonance,  |Y^ 
to 


Re(Yi2}  ^~To  "  +  gC>  (6"44> 

|Re(Y^2)|  and  Re(YC)  =  7T/ncAc,  and  (6-37)  reduces 


T 

res 


4tt 


w  ->  0  ,  ,_a  ,  „c.  2 

wn  n  a  (g  +  g  ) 

a  c  c 


(6-45) 


If  we  introduce  the  parameter 


and  substitute  f or  G  and 
written  as 


v  = 


n  a 

a  a 

n  a 

c  c 


G  as  obtained  from  (6-24), 


T  - 
res  w 


4A  v 

a 


0  w  tt(1  +  v)2 


(6-46) 

then  (6-45)  can  be 

(6-47) 


If  medium  a  and  medium  c  are  the  same,  v  =  1  and  (6-47)  reduces  to 
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(Tw). 


(6-48) 


res  w  -*■  0  it  a 

The  quantity  Tw  is  the  transmission  width,  or  apparent  width,  of  the  slot. 
The  transmission  width  times  the  incident  power  density  equals  the  power 
transmitted  by  a  unit  length  of  the  slot.  Hence,  when  medium  a  and  medium 
c  are  the  same,  at  resonance  the  transmission  width  of  a  narrow  slot  is  1/tt 
wavelengths,  regardless  of  its  actual  width.  When  medium  a  and  medium  c  are 
different,  the  transmission  width  of  a  narrow  slow  at  resonance  can  be  ob¬ 
tained  from  (6-47),  where  v  is  given  by  (6-46).  Note  also  that  (6-47)  and 
(6-48)  are  independent  of  the  medium  in  region  b,  that  is,  the  peak  values 
of  T  are  independent  of  the  medium  in  the  slot.  However,  the  positions 
and  widths  of  the  peaks  are  a  function  of  medium  b,  as  is  evident  from 
(6-43)  with  Y  =  1/wrt. 


6-4  NUMERICAL  RESULTS 

Once  the  transfer  admittance  Y^  is  found  for  a  given  narrow  slot, 
we  can  then  readily  compute  the  field  transmitted  into  the  half-space 
region  c.  It  is  thus  desirable  to  compare  the  results  for  Y^  obtained 
from  the  fequivalent  circuit  formula  (6-29)  with  results  for  Y^  obtained 
from  the  higher  order  moment  solution  of  [4].  In  the  latter,  six  pulse 
expansion  functions  were  used  on  each  aperture  face  and  for  all  slot 
widths  considered.  These  results  are  indicated  at  various  points  on  the 
curves  in  Figs.  6-4  to  6-8  by  markers  (circles,  triangles,  pluses,  and 
crosses).  The  solid  lines  represent  results  obtained  directly  from  the 
equivalent  circuit  model  and  were  computed  a  intervals  of  0.01  wave¬ 
lengths  in  region b.  For  narrow  slots,  the  parameters  Y^  and  T»  computed 
from  the  approximate  solution,  are  periodic  functions  of  d  (screen  thickness). 
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They  are  only  plotted  for  0  <  d/X,  <  2,  where  X.  is  the  wavelength  in 

b  D 

the  slot  region  b.  The  permeability  of  all  regions  is  that  of  free 
space  and  the  permittivity  of  the  different  regions  is  specified  for 
the  different  cases. 

The  real  and  imaginary  parts  of  Y,  „/Y  ,  where  Y  =  1/wri,  is  the 

Iz  o  o  b 

characteristic  admittance  of  the  parallel-plate  transmission  line,  are 

shown  in  Fig.  6-4  for  slot  widths  of  0.01X  to  0.2X  ,  where  X  is  the 

o  o  o 

wavelength  of  free  space.  Equation  (6-28)  was  used  for  the  higher  order 
moment  solution  in  which  the  voltage  becomes  an  integral  over  the 
magnetic  current  M^.  The  two  results  are  in  excellent  agreement  as  w 
becomes  small.  Also,  as  expected  from  (6-29)  and  (6-24),  Y  ^  becomes 
purely  imaginary  as  w  -*  0.  These  computations  indicate  that  the  equi¬ 
valent  circuit  gives  highly  accurate  results  for  w  <  O.IX^,  and  good 

results  for  0.1X  <  w  <  0.2X  .  For  w  >  0.2X  .  the  results  obtained  from 

o  o  o 

the  equivalent  circuit  become  increasingly  inaccurate.  This  is  to  be 
expected,  since  the  half  space  is  then  no  longer  adequately  repre¬ 
sented  by  a  lumped  admittance. 

The  transmission  coefficients  for  the  same  slots  are  shown  in 

Fig.  6-5,  where  the  equivalent  circuit  results  were  obtained  from  (6-37). 

For  the  higher  order  moment  solution,  P  was  obtained  by  the  usual 

trans 

multi-port  network  formula  [4] .  As  w  -*  0,  T  becomes  a  maximum  at  so- 

called  resonant  thicknesses,  which  approach  multiples  of  X. / 2 ,  as  pre- 

fa 

dieted  by  (6-42).  The  peak  value  of  T  for  w  =  0.0lXo  in  Fig.  6-5  is 
T  =  30.5  occurring  at  d/Xfe  =  0.47,  0.97,  1.47,...  .  This  is  in  agree¬ 
ment  with  the  prediction  that  the  transmission  width  Tw -*  X  /it  at  values  of  d 

o 

slightly  less  than  integral  multiples  of  X,  /I  as  w  0. 

b 
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Fig.  6-5.  Plots  of  transmission  coefficient  T  vs.  d/A  for  cases  of 

b 

Fig.  6-4.  Width  w  =  0.01A  (circles),  w  =  0.05A  (triangles), 

o  o 

w  =  0.1A  (pluses),  and  w  =  0.2A  (crosses), 
o  o 

When  the  slot  is  loaded  with  different  dielectrics,  the 

imaginary  part  of  Y  ^  is  most  affected.  This  is  shown  in  Fig.  6-6  for  a 

slot  of  width  w  =  0.1A  for  e,  =  e  ,  5e  ,  and  10e  ,  where  c  is  the 

o  b  o  o  o  o 

permittivity  of  free  space.  As  the  material  filling  the  slot  becomes  more 

dense,  the  transmission  resonances  become  narrower  and  occur  closer  to 

multiples  of  A,  /2.  This  is  shown  in  Fig.  6-7  for  the  slot  of  Fig.  6-6, 
b 

the  result  being  in  agreement  with  (6-43).  Note  that  the  heights  of  the 
peaks,  however,  are  the  same. 


Fig.  6-6.  Plots  of  Re(Y.„/Y  )  and  Im(Y../Y  )  vs.  d/\v  for  w  =  0.1 

■L^  o  1/  o  b 

ea  m  ec  m  e0»  wlth  various  dielectrics  filling  the  slot 
Cases  shown  are  efe  -  zq  (circles),  eb  =  5eQ  (triangles) 
and  e.  ■  10A  (pluses) . 
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* 


Fig.  6-7.  Plots  of  transmission  coefficient  T  vs.  d/A^  for  the 

case-"  of  Fig.  6-6.  The  dielectrics  filling  the  slot 

are  e,  =  e  (circles),  e.  =  5e  (triangles),  and 
bo  bo 

e,  =  lOe  (pluses). 

D  O 

The  effects  of  lossy  material  filling  the  slot  are  shown  in 

Fig.  6-8.  Here  the  slot  width  is  0.05Ao.  The  expected  decay  in  the 

transmission  coefficient  peaks  is  seen  to  occur  with  increasing  screen 

thickness  d.  Circles  represent  the  lossless  case  e.  =  c  ,  triangles 

b  o 

represent  a  dielectric  e,  =  (l-jO.Ol)t  ,  that  is,  with  Q  factor  =  100, 

b  o 

and  pluses  represent  a  dielectric  =  (l-j0.1)£o,  that  is,  with  Q  factor 
=  10.  Note  that  loss  does  not  affect  the  position  of  the  resonances. 

If  the  conductor  has  large  but  finite  conductivity,  the  effect  of 
the  conductor  loss  is  similar  to  that  of  dielectric  loss.  For  small 


Fig.  6-8.  Plots  of  transmission  coefficient  T  vs.  d/A,  for  f  =  c  =  f  , 

b  a  c  o 

w  =  0.05Ao,  for  lossy  dielectric  in  the  slot.  Cases  shown  are 

=  eq  (circles),  =  (l-j0.01)eo  (triangles),  and 

e,  =  (l-jO.l)e  (pluses), 
b  o 

losses,  either  conductor  or  dielectric,  the  principal  effect  is  to  change 
the  propagation  constant  of  the  parallel-plate  transmission  line  mode  from 
purely  imaginary  to  complex  with  a  ?•  ’ '  eal  part  [11].  Equating  the 

attenuation  constant  for  a  parallel-plate  iveguide  with  perfectly  con¬ 
ducting  walls  and  lossy  dielectric  (e  =  e'  -  je"  and  quality  factor 
Q  *  e'/c")  to  one  with  perfect  dielectric  and  lossy  conducting  walls 

(surface  resistance  R  =  vwy/2 a 7),  we  obtain 

s  cond 
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Now  the  lossy  conductor  case  behaves  similarly  to  the  lossy  dielectric 
case  with  Q  factor  given  by  (6-49) .  Note  that  Q  is  proportional  to  w, 
meaning  that  for  very  narrow  slots  the  losses  due  to  conducting  walls 
become  an  important  factor  limiting  the  field  penetration  through  slots 
in  conductors  of  resonant  thickness. 

6-5  CONCLUSION 

The  validity  of  an  equivalent  circuit  model  for  a  narrow  TE  excited 
slot  in  a  thick  conducting  screen  has  been  investigated.  The  equivalent 
circuit  formulas  for  transfer  admittance  and  transmission  coefficient  give 
accurate  results  when  compared  to  a  higher  order  moment  solution.  When 
the  material  filling  the  slot  is  dense  enough  to  allow  more  than  one  propa¬ 
gating  mode,  the  equivalent  circuit  picture  of  Fig.  6-3  becomes  more  com¬ 
plicated.  One  must  then  consider  a  sequence  of  transmission  lines  con¬ 
necting  regions  a  and  c  in  which  the  higher  order  modes  couple  to  one 
another . 

An  interesting  result  of  this  investigation  is  that  Tw,  the  trans¬ 
mission  width  of  the  slot,  for  a  narrow  slot  at  resonance,  is  independent 
of  actual  slot  width.  This  phenomenon  is  analogous  to  the  phenomenon  of 
scattering  by  resonant  scatters,  or  reception  by  resonant  antennas.  For 
example,  a  short  dipole  resonated  by  an  inductor  has  a  scattering 
cross  section  of  9/ Att  square  wavelengths  regardless  of  the  actual  size  of 
the  dipole  [  12  ] .  A  short  dipole  resonated  by  an  inductor  and  used  as  a 
receiving  antenna  has  an  effective  aperture  of  3/8tt  square  wavelengths 
regardless  of  the  actual  size  of  the  dipole  f 13 ] . 
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Chapter  VII 

RESONANT  BEHAVIOR  OF  A  SMALL  APERTURE 
BACKED  BY  A  CONDUCTING  BODY 

7-1  INTRODUCTION 

The  coupling  of  electromagnetic  energy  to  wires  or  other  conducting 
objects  through  an  aperture  in  a  conducting  wall  is  an  important  problem 
in  the  theory  of  electromagnetic  compatibility  and  interference.  As  an 
approximation,  some  investigators  have  first  solved  the  problem  for  the 
tangential  electric  field  in  the  aperture,  and  then  taken  this  as  a 
secondary  source  for  the  field  coupled  to  a  wire  [1] .  This  approach 
neglects  the  effect  of  the  wire  on  the  electric  field  in  the  aperture, 
which  can  be  appreciable.  Under  certain  conditions  of  resonance,  the 
power  transmitted  through  an  aperture  can  be  orders  of  magnitude  larger 
when  an  object  is  near  it  than  when  no  object  is  present. 

Our  approach  is  to  first  obtain  the  functional  equations  for  the 
problem  using  the  equivalence  principle  (Sec.  3-5  of  [2]),  and  then  to 
reduce  these  equations  to  matrix  form  via  the  method  of  moments  [3] .  The 
various  matrices  are  interpreted  in  terms  of  generalized  network  parame¬ 
ters,  such  as  voltages,  currents,  admittances,  and  impedances  [3].  The 
aperture  admittance  matrices  of  electrically  small  apertures  are  obtained 
and  discussed.  An  example  of  coupling  through  a  capacitively  loaded 
aperture  is  given  to  illustrate  the  phenomenon  of  aperture  resonance. 

Perhaps  the  first  accurate  treatment  of  coupling  through  an  aper¬ 
ture  to  a  wire  was  given  by  Butler  and  Umashankar  [4] .  Later  work  by 
Butler  considered  the  electrically  small  aperture  and  coupling  to 
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objects  other  than  wires  [5].  He  used  Bethe-hole  theory  for  the  small 
aperture  but  did  not  include  a  radiation  term.  He  also  did  not  look  at 
the  resonance  effects  which  we  emphasize.  Other  work  which  treats 
coupling  to  wires  is  that  of  Kajfez  [6],  Kajfez  and  Wilton  [7],  and  lee 
and  Yang  [8].  These  last  three  references  deal  with  determining  an 
equivalent  circuit  for  wires  passing  near  small  apertures.  Again  Bethe- 
hole  theory  was  used,  radiation  was  not  accounted  for  in  the  equivalent 
circuit,  and  resonance  effects  were  not  considered. 

7-2  FORMULATION  OF  THE  PROBLEM 

The  general  problem  of  coupling  to  a  conducting  body  through  an 
aperture  in  a  conducting  wall  is  represented  by  Fig.  7-1.  The  wall 


CONDUCTING 

BODY 


Fig.  7-1. 


Two-region  aperture  problem  with  impressed  sources  J*,m' 
in  region  a  and  conducting  object  in  region  b. 
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divides  space  into  two  regions,  called  region  a  and  region  b.  The  exci¬ 
tation  is  represented  by  impressed  sources  in  region  a,  and  the  conduct¬ 
ing  body  is  in  region  b.  The  problem  is  primarily  that  of  finding  the 
tangential  electric  field  in  the  aperture  and  the  current  on  the  con¬ 
ducting  body,  and  secondarily  that  of  finding  the  fields  throughout  space. 

We  use  the  equivalence  principle  (Section  1-6)  to  divide  the  problem 
into  two  parts,  as  shown  in  Fig.  7-2.  The  field  in  region  a  remains 
unchanged  if  the  aperture  is  closed  by  a  conductor  and  the  equivalent  sur¬ 
face  magnetic  current 

M  =  n  x  E  (7-1) 

is  placed  over  the  aperture  region,  where  E  is  the  electric  field  in  the 
aperture  of  the  original  problem  and  n  is  the  unit  normal  pointing  out¬ 
ward.  The  original  sources  J*,  M*'  must  be  kept  in  region  a.  This  equiva¬ 
lence  is  shown  in  Fig.  7-2a.  The  field  in  region  b  remains  unchanged  if 
the  aperture  is  closed  by  a  conductor  and  the  equivalent  surface  magnetic 
current  -M  is  placed  over  the  aperture  region.  The  current  J  on  the  con¬ 
ducting  body  must  be  kept  in  region  b.  This  equivalence  is  shown  in 
Fig.  7-2b. 

The  fact  that  M  is  used  in  Fig.  7-2a  and  -M  is  used  in  Fig.  7-2b 
ensures  that  the  tangential  electric  field  is  equal  on  each  side  of  the 
aperture  region.  We  have  two  more  boundary  conditions  to  enforce: 

(a)  the  tangential  magnetic  field  must  be  equal  on  each  side  of  the  aper¬ 
ture  region,  and  (b)  the  tangential  electric  field  must  vanish  on  the 
surface  of  the  conducting  body.  These  two  conditions  give  us  two  equa¬ 
tions  from  which  to  calculate  the  unknown  quantities  M  and  J. 


I 


REGION  a 


f  CONDUCTOR 


(a)  EQUIVALENCE  FOR  REGION  o 


CONDUCTOR 
(APERTURE  SHORTED) 


lb)  EQUIVALENCE  FOR  REGION  b 


Fig.  7-2.  Original  problem  divided  into  two  equivalent  parts. 


To  express  these  boundary  conditions  in  equation  form,  we  intro¬ 
duce  the  following  notation:  Let  the  electric  and  magnetic  fields  in 
region  a  be  denoted 

E3  =  Ea(M)  +  E1  (7-2 

H3  =  Ha(M)  +  H1  (7-3 

where  Ea(M),  Ha(M)  are  the  fields  from  M  in  Fig.  7-2a,  and  E*,  are 
the  fields  from  J^-,  in  Fig.  7-2a.  Note  that  all  fields  are  computed 
with  the  aperture  shorted.  Similarly,  let  the  electric  and  magnetic 
fields  in  region  b  be  denoted 

Eb  =  Eb(-M)  +  Eb(J)  (7-4 

Hb  =  Hb(-M)  4-  Hb(J)  (7-5) 

K  K 

where  j£  (-M) ,  H  (-M)  are  the  fields  from  -M  in  Fig.  7-2b,  and  ET(J), 

Hb(J)  are  the  fields  from  J_  in  Fig.  7-2b.  Again,  all  fields  are 
computed  with  the  aperture  shorted. 

To  satisfy  the  boundary  condition  that  the  tangential  component 
of  H  must  be  continuous  across  the  aperture,  we  equate  tangential  com¬ 
ponents  of  (7-3)  and  (7-5),  obtaining 

-HJJ(M)  -  H^(M)  +  hJ(J)  -  +  hJ  over  A  (7-6) 

where  A  denotes  the  aperture  region.  The  subscripts  t  denote  the  tan¬ 
gential  component  over  A,  and  we  have  used  the  linearity  of  the  operator 
to  replace  H^(-M)  by  -H^(M).  To  satisfy  the  boundary  condition  that  the 
tangential  compoi  ent  of  E  must  vanish  on  the  surface  of  the  conducting 
body,  c -*t  i  tangential  component  of  (7-4)  equal  to  zero,  obtaining 
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Eb(M)  -  E^(J)  =0 


over  B 


(7-7) 


where  B  denotes  the  conducting  body  surface.  The  subscripts  t  denote 
the  tangential  component  over  B,  and  we  have  used  the  linearity  of  the 
operator  to  replace  E^(-M)  by  -E^ (M) .  Equations  (7-6)  and  (7-7)  are 
vector  equations  for  determining  the  unknowns  M  and  which  exist  over 
the  surfaces  A  and  B. 

We  next  reduce  (7-6)  and  (7-7)  to  matrix  equations  using  the  method 

of  moments  [3].  For  this,  we  define  a  set  of  expansion  functions 

(M  ,  n  =  1,2,...,  N.  1  and  express  the  magnetic  current  over  A  as 
— n  A 


M  =  y  V  M  (7-81 

—  u  n~n 
n 


where  V  are  coefficients  to  be  determined.  We  define  a  set  of  expan- 
n 

sion  functions  {j  ,  n  -  1,2,.,.,  N  }  and  express  the  electric  current 
— n  B 

over  B  as 


J 


(7-9) 


where  I  are  coefficients  to  be  determined.  We  substitute  (7-8)  and 
n 

(7-9)  into  (7-6)  and  (7-7),  use  the  linearity  of  the  operators,  and 
obtain 


(7-10) 


over  A,  and 


Y  v  eNm  )  -  l  I  eVi  )  =  0 

L  n— t  — n  L  n~t  — n 

n  n 


(7-11) 


over  B.  For  A,  we  define  a  set  of  testing  functions  ,  m=l,2,. 

-~m 


and  a  symmetric  product 


<F,  G'>  = 

A  ) 


F  •  G  ds 


(7-121 
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We  take  the  symmetric  product  of  (7-10)  with  each  M  ,  and  use  the  linearity 

- m 

of  the  symmetric  product  to  obtain 


y  v  <ft  , 

u  n  m 
n 


H*(M  )> 
t  — n  A 


y  v 

*-»  n  tt 


(7-13) 


+  y  I  <M  ,  Hb(j  )>.  =  ,  H1>. 

°  n  m  tv  n  A  m  t  A 
n 


m  =  1,2,...,  N^.  For  B,  we  define  a  set  of  testing  functions 

{j  ,  m  =  1,2,...,  N_}  and  a  symmetric  product 
— m  d 


<F, 


jj  F  •  G  ds 
B 


(7-14) 


We  take  the  symmetric  product  of  (7-11)  with  each  J^,  and  use  the  linearity 
of  the  symmetric  product  to  obtain 


l  V  <J  ,  Eb(M  )>  -  l  I  <J  ,  Eb(J  )>  *  0 

L  n  m  t  n  B  Ln  m  t  n  B 


(7-15) 


m  =  1,2,...,  N  .  Equations  (7-13)  and  (7-15)  are  now  a  set  of  algebraic 

D 

equations  for  determining  the  unknown  coefficients  and  1^. 

The  above  equations  can  be  put  into  matrix  notation  as  follows: 
Define  an  admittance  matrix  for  region  a  as 


A  A 


and  an  admittance  matrix  for  region  b  as 


[vh  -  [<-«,, 

A  A 


(7-16) 


(7-17) 


Note  that  these  are  exactly  the  same  as  defined  in  Chapter  II  for  the 
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problem  with  the  wire  absent.  Define  coupling  matrices 


[T!  -  [<ftn,  Ht<Jn)>A1  „  „ 

A  D 


[T]  -  [<Jm,  Et<Mn)>BlN  VM 

B  A 

Define  an  impedance  matrix  for  the  wire  object  as 


iz-  ■  tH,-  WVn  *». 

D  O 

Note  that  this  is  calculated  in  the  presence  of  a  complete  conducting 
boundary,  that  is,  with  the  aperture  shorted.  Define  a  source  vector 

11  ‘  [<V  "^aV  x  1 

A 


(7-18) 


(7-19) 


(7-20) 


(7-21) 


and  coefficient  vectors 


V  '  tVn'»A  «  1 

A 


1  ■ 


(7-22) 


(7-23) 


Now  the  matrix  equations  equivalent  to  (7-13)  and  (7-15)  are 


[Ya  +  Yb]V  +  [T]I  =  I1 


(7-24) 


[T]  V  +  [Z]I  =  0  (7-25) 

The  coefficient  vectors  V  and  I  are  obtained  from  the  matrix  solution 
to  (7-24)  and  (7-25).  If  we  want  only  I,  we  can  eliminate  V  from  (7-25) 
by  using  (7-24).  Similarly,  if  we  want  only  V,  we  can  eliminate  I  from 
(7-24)  by  using  (7-25). 
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For  emphasis,  we  restate  that  [Y*3]  is  calculated  for  the  aper¬ 
ture  as  if  the  conducting  body  were  not  present,  by  methods  outlined  in 
Chapter  II  and  carried  out  in  detail  for  a  rectangular  aperture  in  a 
plane  conducting  wall  in  Chapter  III.  The  matrix  [Z]  is  calculated  for 
the  conducting  body  as  if  the  aperture  were  not  present,  tiiat  is,  with 
the  aperture  short  circuited.  Conceptually,  the  only  new  matrices  re¬ 
quired  are  the  interaction  matrices,  [T]  and  [T].  If  a  Galerkin  solu¬ 
tion  is  used,  that  is,  if  (M  }  =  {M  }  and  {j  }  =  {J  },  it  then  follows 

~n  -n  — n  — n 

from  reciprocity  (Sec.  3-8  of  [2])  that  [T]  =  -  [T],  where  the  tilde 
denotes  transpose. 

7-3  APERTURE  ADMITTANCE  MATRIX  FOR  SMALL  HOLES 

An  aperture  admittance  matrix  for  electrically  small  apertures 
can  be  obtained  from  the  quasi-static  solution  to  the  integral  equation. 
The  complete  solution  depends  on  both  regions  a  and  b,  but  the  components 
of  the  aperture  admittance  matrix  depend  only  on  one  region.  The  usual 
procedure  is  to  solve  for  the  aperture  admittance  of  a  canonical  problem, 
and  then  use  this  solution  for  other  problems  which  differ  only  slightly 
from  the  canonical  problem.  The  basic  canonical  problem  for  aperture 
admittance  is  that  of  an  aperture  in  an  infinite  plane  conducting  screen, 
excited  by  an  incident  plane  wave.  The  electrostatic  and  magnetostatic 
treatment  of  this  problem  is  known  as  Bethe-hole  theory  [10],  [11]. 

The  field  produced  by  an  aperture  is  defined  to  be  the  dif¬ 
ference  between  that  produced  by  the  original  sources  in  the  presence 
of  a  shorted  aperture  and  that  produced  when  the  aperture  is  present. 
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Except  in  the  immediate  vicinity  of  the  aperture,  the  field  of  an  elec¬ 
trically  small  aperture  in  a  conducting  screen  is  that  of  a  magnetic 
dipole  tangential  to  the  screen  plus  that  of  an  electric  dipole  £ 
normal  to  the  screen,  both  radiating  with  the  aperture  shorted.  Let 
the  two  half  space  regions  on  each  side  of  the  screen  be  denoted  region 

S  C  cl  S  C  cl 

a  and  region  b,  with  sources  possibly  in  both  regions.  Let  H  ,  E, 
set)  seb 

and  H  ,  E!  be  the  fields  produced  by  these  sources  when  the  aperture 
is  shorted.  Then  the  difference  field  in  region  b  is  given  by  that  from 
the  dipoles 

=  ~  %  *  (E  ~  H  )  (7-26) 

2^  =  e  ag  •  (E  -  E  )  (7-2/ ) 

where  am  and  are  the  tensor  magnetic  and  electric  polarizabilities 


a  =  a  ,  t,  t.  +  a  -,t0t0 
m  ml— 1—1  m2— 2— 2 


(7-28) 


a  =  a  n  n  (7-29) 

e  e  —  — 

Here  t^  and  _t  are  unit  vectors  tangential  to  the  screen  chosen  to 

diagonalize  a^,  and  ri  is  the  unit  normal  to  the  screen  pointing  into 

region  b.  Analytical  solutions  for  m  and  a  exist  for  circular  aper- 

m  e 

tures  [10]  and  elliptical  apertures  [12].  Methods  for  computing  <m  and 

for  other  shapes  have  been  developed  by  De  Meulenaere  and  Van  Blade]  [13], 
and  of  measuring  them  by  Cohn  [14],  [15].  For  an  aperture  with  an  axis  of 
symmetry,  _t^  is  in  the  direction  of  this  axis  and  t_2  is  normal  to  it. 


106 


The  magnetic  dipole  moment  is  that  of  a  magnetic  pole  dipole 
and  the  electric  dipole  moment  is  that  of  an  electric  charge  dipole. 

For  our  purposes  it  is  more  convenient  to  deal  with  current  elements, 
magnetic  denoted  by  K2.  and  electric  denoted  by  I_£.  These  are  related 
to  and  by 

=  jwy  2^  (7-30) 


U  =  jw 


(7-31) 


scci  sea  i 

Finally,  let  sources  exist  only  in  region  a,  and  denote  j2  . ,  H  by  E  , 
H1,  as  was  done  in  section  7-2  or  in  Chapter  II.  Then  (7-30)  and  (7-31), 
with  substitutions  from  (7-26)  and  (7-27),  become 


K£  =  -  jcoy  am  ♦  H  (0) 


(7-32) 


U  =  jwe  ag  •  E  (0) 


(7-33) 


i  i  i  i 

The  notation  E_  (0),  H  (0)  denotes  the  value  of  12  and  H  at  the  center  0 

of  the  aperture.  These  impressed  fields  are  assumed  to  be  approximately 

constant  over  the  aperture  region. 

We  wish  to  translate  the  above  results  into  those  of  aperture 

admittance,  and  to  augment  them  to  include  the  effects  of  radiation. 

Hence,  for  a  small  aperture,  let  the  magnetic  current  be  expressed  by  .a 

three  term  expansion 


M  =  V1M1  +  V2M2  +  V3M3 


(7-34) 


where  is  the  quasi-static  current  which  produces  the  effect  of  a 
unit  magnetic  dipole  K£  =  1  in  the  direction,  *s  t*ie  quasi-static 


E 


current  which  produces  the  effect  of  a  unit  magnetic  dipole  K2.  =  1  in 
the  t^  direction,  and  is  the  quasi-static  current  which  produces  the 
effect  of  an  electric  dipole  1)1  3  j  in  the  n  direction.  (A  real  mag¬ 
netic  current  gives  rise  to  an  imaginary  12.  [2,  p.  135].)  We  will  not 
need  explicit  formulas  for  M^,  M2,  and  M^,  but  they  can  be  obtained  from 
the  Bethe-hole  quasi-static  solution. 

The  dipoles  in  (7-32)  and  (7-33)  radiate  in  region  b.  In  Fig.  7-2b, 
-M  radiates  in  region  b.  The  dipole  effects  of  -M  will  be  equal  to  the 
results  (7-32)  and  (7-33)  of  Bethe-hole  theory  if  the  V's  in  (7-34)  are 
given  by 

vi  - 

V2  =  ju)uam2H2(0)  (7-35) 

V,  =  -  wea  E*(0) 

3  e  i 

Here,  the  subscripts  1,  2,  and  3  denote  the  t^,  t_2>  and  n  directions, 
respectively.  For  our  canonical  problem  involving  the  infinite  plane 
conducting  screen,  1  m  <5  in  (7-24)  so  that  the  method  of  moments  solution 
for  ^  reduces  to 

V  -  [Ya  +  yh]'1  I1  (7-36) 

For  testing  functions  we  use  n  =  1,2,3  (Galerkin's  method).  We 


then  find 
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Now  for  the  admittance  solution  (7-36)  to  give  the  same  result  as  Bethe 
hole  theory,  we  must  have 


[Y  +  Y  ]  = 


jwyct 


ml 


0 


0 


1 

jWuam2 

0 


0 

-1 

jweae 


(7-38) 


The  solution  (7-36)  gives  coefficients  which,  according  to  (7-34),  give 
the  negative  of  the  dipoles  specified  by  (7-32)  and  (7-33). 

The  aperture  admittance  matrix  obtained  from  Bethe-hole  theory  neg¬ 
lects  radiation.  In  other  words,  it  evaluates  only  the  first  term  of  a 

frequency  expansion  for  [Y3  +  Yb] .  When  there  is  an  object  near  the  aper- 

3  b 

ture,  as  in  our  original  problem,  the  first  (susceptive)  term  in  [Y  +  Y  ] 
can  be  cancelled  by  the  interaction  between  the  object  and  the  aperture, 
resulting  in  the  incorrect  prediction  of  infinite  coefficients  V^,  or  infi¬ 
nite  power  transmitted  through  the  aperture.  This  defect  in  the  Bethe- 

hole  theory  can  be  corrected  by  evaluating  the  second  (conductive)  term  of 

g  b 

a  frequency  expansion  for  [Y  +  Y  ].  Fortunately,  this  additional  term  is 

easily  obtained  from  the  radiation  field  of  the  dipoles. 

3  b 

Note  that,  for  M  real,  we  have  Re(Y  +  Y  )  equal  to  4  times 
— n  nn  nn 

the  power  radiated  by  in  free  space.  The  factor  of  4  arises  because 

H^M^)  in  (7-16)  is  twice  the  magnetic  field  due  to  radiating  in  free 

space  and  because  [Ya  +  Yb]  =  2[Ya].  The  field  distant  from  M  is  assumed 

— n 

to  be  a  dipole  field,  hence  we  can  use  the  results  for  a  dipole.  Using 
duality  (Sec.  3-2  of  [2])  and  the  formula  for  the  power  radiated  by  an 
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electric  dipole  (Eq.  2-116  of  [2]),  we  have  for  a  magnetic  dipole 

1 2 


P  =  ^ 
3n 


K£ 


(7-39) 


where  n  *  /y/e  is  the  intrinsic  impedance  of  space  and  X  is  the  wave¬ 
length.  Since  and  M2  each  produce  the  effect  of  Ki  »  1,  we  have 


Re[Ya  +  Yb]u  -  Re[Ya  +  Yb]22  =  4P 


8tt 

3nX2 


(7-40) 


By  similar  reasoning,  for  the  electric  dipole  l£  =  j  radiating  on  the 
surface  of  an  infinite  conducting  plane  we  obtain  a  term 


Re[Ya  +  Yb]  =  ^ 
3XZ 


(7-41) 


These  terms  must  be  added  to  the  imaginary  terms  in  (7-38)  to  give  the 
aperture  admittances  corrected  to  account  for  radiation.  The  off-diagonal 
elements  in  (7-38)  remain  zero  because  there  is  no  cross-power  between  any 
pair  of  the  quasi-static  currents  M^,  M2>  and  M^. 

For  a  half-space  region  alone  the  matrix  elements  are  one-half  those 
evaluated  above,  or 


[Y 


hs. 


where 


hs 

11 

0 

0 

0 

hs 

22 

0 

0 

0 

hs 

33 

hs 

1 

X 

4tt 

11  = 

2ju>yami 

T 

3nX2 

hs 

1 

1 

4tt 

22 

2ju>li«m2 

▼ 

2 

3nx 

hs 

-1 

X 

4TTH 

33  “ 

2  ju)Cae 

• 

3X2  . 

(7-42) 


(7-43) 
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Of  course,  equations  (7-43)  are  strictly  valid  only  for  a  half  space 
region,  but  they  are  often  a  good  approximation  for  other  regions. 


7-4  TRANSMITTED  POWER 


The  complex  power  through  the  aperture  is  basically  given  by 


P 


A 


j  j  E  x  H*  •  n  ds 
A 


jj  M  •  H*  ds 
A 


(7-44) 


where  M  is  the  equivalent  magnetic  current  defined  by  (7-1),  and  H*  is 
the  conjugate  of  the  magnetic  field  in  the  aperture.  Setting  H  =  Hb, 
and  using  (7-5) ,  we  have 


r  r 


[-Hb(M)  +  Hb(J)1*  ds 


b* 

H  (M)ds 


M 


Hb*( J) 


ds 


(7-45) 


Substituting  for  M  from  (7-8)  and  for  .J  from  (7-9),  we  reduce  this  to 


P.  =  "  I  l  V  V* 
A  u  u  m  n 
m  n 


+  n  v  i* 

L  u  m  n 

m  n 


M  ♦  Hb*(M  )  ds 
-m  —  — n 


M  •  Hb*(J  )  ds 


(7-46) 


If  a  Galerkin  solution  (M  =  M  and  J  =  J  )  is  used,  and  if  the  expansion 

— n  — n  -n  — n 
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functions  are  real,  then  the  negatives  of  the  first  integrals  in  (7-46) 

b* 

are  Y  as  defined  by  (7-17),  and  the  second  integrals  in  (7-46)  are  T* 
mn  mn 

as  defined  by  (7-18).  Hence, 


l  l 


b* 

V  V*Y 
m  n  mn 


I  I  V  I*T* 

L  L  m  n  mn 


In  matrix  form  this  becomes 


(7-47) 


PA  *  V f Yb ] *V*  +  V[T]*I*  (7-48) 

where  V  is  defined  by  (7-22),  I  by  (7-23),  and  the  tilde  denotes  transpose. 

The  last  term  in  (7-48)  can  be  put  into  another  form  by  using  (7-25) 
to  express  t  in  terms  of  Since  [f]  “  -  [T],  we  have 

I  =■  [Z_1T]V  (7-49) 

Substituting  this  into  (7-48),  and  combining  the  two  matrix  products,  we 
have 

PA  ”  ^[Yb  +  TZ~lTl*  V*  (7-50) 

This  suggests  that  we  can  define  the  effective  aperture  admittance  into 
region  b  to  be 

[Yeff!  ’  [Yb  +  TZ~1yI  (7_51) 

Then  (7-50)  becomes  the  usual  formula  for  power  transmitted  into  an  N-port 
network.  Note  that  the  effective  aperture  admittance  is  that  for  region  b 
with  the  conducting  body  present.  In  this  case  we  are  viewing  the  con¬ 
ducting  body  as  part  of  the  definition  of  region  b,  that  is,  part  of  the 
environment  into  which  the  aperture  radiates. 


We  can  see  more  clearly  what  can  happen  if  the  magnetic  current 
M  in  the  aperture  and  the  electric  current  J  on  the  body  are  expressed 


in  terms  of  a  single  expansion  function  each.  These  two  quantities 
are  then  given  by 

M  =  VM  and  J  =  IJ^ 

where  only  the  complex  amplitudes  V  and  I  are  unknown.  The  matrix 
equations  (7-24)  and  (7-25)  are  now  reduced  to  the  scalar  equations 

YV  +  TI  *  H®Ca 

-TV  +  ZI  =  0 

where  Y  =  Ya  +  y'5,  T,  and  Z  are  given  by  (7-16,  17,  18,  and  20)  with 
the  subscripts  m  and  n  each  being  1  only.  Expressing  I  in  terms  of 
V  from  (7-54),  and  substituting  into  (7-53),  we  have 

(Y  +  y-)V  =  Haca  I 

For  electrically  small  apertures,  normally  |lm(Y)|  »  Re(Y).  We  define 
aperture-body  resonance  to  be  the  case  for  which 

Im(T2/Z)  -  -  Im(Y) 

2 

T  is  real  positive  if  the  body  and  aperture  are  electrically  small 
and  close  to  each  other.  For  this  case,  at  resonance  (7-55)  reduces  to 

<G  +  — j)  V  -  H®c*  I 

|z|  r  * 

where  the  subscript  r  on  denotes  "at  aperture-body  resonance," 

G  *  Re(Y) ,  and  R  *  Re(Z). 
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In  the  next  section  we  show  by  example  that  the  power  transmitted 
by  an  electrically  small  aperture  backed  by  a  conducting  body  can  be  much 
larger  when  the  body  is  present  than  when  the  body  is  absent.  As  a  limit¬ 
ing  case,  the  second  term  in  parenthesis  in  (7-55)  may  be  small  compared 
to  the  first  term.  When  the  aperture  is  far  from  resonance,  we  then  have 


V 


.sea 


Im(Y) 


(7-58) 


When  aperture-body  resonance  occurs,  we  have 


H 


sea 


Re(Y) 


We  define  the  aperture  Q  as 


llm(Y) 

Re(Y) 


(7-59) 


(7-60) 


Then,  since  M  *  VM, ,  we  see  that  in  the  aperture  can  be  increased  at 
—  —I  — t 

most  by  the  factor  Q  at  aperture-body  resonance.  The  power  transmitted 

2 

by  the  aperture,  (7-50),  depends  on  |v|  ,  and  hence  may  be  increased  at 
2 

most  by  Q  .  As  a  word  of  caution,  note  that  Q  depends  on  as  well  as  on 
the  aperture  size  and  shape. 


7-5  AN  EXAMPLE 

As  an  example,  consider  a  small  rectangular  aperture  in  a  plane 
conducting  screen  with  a  capacitor  placed  across  its  midpoint,  as  shown 
in  Fig.  7-3a.  Take  the  excitation  to  be  a  plane  wave  normally  incident 
on  the  screen  with  E*  perpendicular  to  the  longer  axis  of  the  slot.  By 
Babinet's  principle  (Sec.  7-12  of  [2]),  this  problem  is  dual  to  that  of 
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Fig.  7-3.  (a)  A  capacitively  loaded  aperture  and  (b)  the 

complementary  inductively  loaded  scatterer. 


electromagnetic  scattering  from  a  conducting  rectangular  dipole,  with 
an  inductor  in  series  at  its  midpoint,  as  shown  in  Fig.  7-3b.  It  is 
known  that,  for  a  loss-free  resonated  dipole,  the  back  scattering  cross 
section  is 


a 


(7-61) 


regardless  of  the  size  or  shape  of  the  dipole  [16].  The  result  (7-61) 
can,  however,  be  greatly  decreased  by  conductor  losses  for  small  dipoles. 
We  now  show  that  a  similar  result  holds  for  the  transmission  cross  section 
of  a  small  resonated  aperture. 


For  a  plane  wave  normally  incident  on  an  aperture,  we  define  the 
transmission  coefficient  T  of  the  aperture  to  be  the  ratio  of  the  power 
transmitted  through  the  aperture  to  the  power  incident  on  it,  that  is. 


115 


(7-66) 
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SCH 

Finally,  using  (7-64)  and  the  fact  that  =  2Hq,  we  obtain 


T  = 


hag 


_  31 
b  47TA 


(7-67) 


Hence,  the  transmission  coefficient  of  an  electrically  small  aperture 
resonated  by  a  capacitor  is  independent  of  the  shape  of  the  aperture. 

The  transmission  area  of  the  aperture  is  defined  as  the  area  for  which 
the  incident  wave  contains  the  same  power  as  transmitted  by  the  aperture, 
that  is 


T  A  = 


3A 

4ir 


(7-68) 


Hence,  the  transmission  area  of  a  small  aperture  resonated  by  a  capacitor 
2 

is  3A  /4tt  independent  of  the  size  and  shape  of  the  aperture.  Of  course, 
conductor  losses  may  greatly  decrease  (7-68)  for  very  small  apertures. 


7-6  DISCUSSION 

This  chapter  illustrates  how  a  conducting  object  situated  near  a 

small  aperture  in  a  conducting  screen  can  resonate  the  aperture,  thereby 

increasing  the  power  transmitted  by  that  aperture  over  what  it  would 

transmit  with  no  object  present.  The  magnitude  of  the  increase  can  be  of 
2 

the  order  of  Q  ,  where  Q  is  the  quality  factor  of  the  aperture  current 
being  resonated.  With  no  conductor  losses  Q  becomes  very  large  for 
electrically  small  holes.  For  example,  the  magnetic  dipole  of  a  circular 
aperture  of  radius  R  has  a  Q  of  114  when  R  =  A/20  and  a  Q  of  14,200  when 
R  =  A/100.  These  values  of  Q  were  obtained  from  (7-60),  (7-43),  and, 
as  given  in  [51,  a  ^  =  (4/3)R  .  When  losses  are  included,  we  can  expect 


Q  to  be  limited  to  the  order  of  a  few  hundred  at  radio  frequencies. 
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The  susceptance  of  the  magnetic  dipole  mode  for  a  small  aperture 
is  inductive,  and  requires  the  coupling  of  a  capacitive  susceptance  from 
the  backing  object  for  resonance.  In  the  example  chosen,  we  obtained 
this  required  susceptance  by  capacitively  loading  the  slot.  We  used  a 
lumped  load,  but  longer  wires  can  produce  this  capacitive  susceptance 
without  lumped  loads.  For  example,  a  straight  wire  of  length  l  slightly 
less  than  a  half  wavelength  (a  wire  dipole)  produces  a  capacitive  sus¬ 
ceptance.  If  region  b  is  a  cavity,  it  will  reflect  a  capacitive  sus¬ 
ceptance  at  some  frequencies  and  an  inductive  susceptance  at  other  fre¬ 
quencies.  The  reason  we  chose  a  lumped  capacitive  load  was  to  show  that 
aperture-conductor  resonance  is  possible  even  if  all  dimensions  are 
electrically  small. 

Another  point  that  we  wish  to  emphasize  is  that  Bethe-hole  theory 
should  assume  that  the  form  of  the  aperture  field  remains  almost  the  same 
as  in  the  canonical  problem  of  a  small  aperture  in  a  plane  ^  teen,  but 
not  its  amplitude  and  phase.  If  this  is  not  recognized,  order  of  magni¬ 
tude  errors  may  result  when  the  aperture  interacts  with  other  objects. 

Also,  the  concept  of  small-aperture  polarizability  should  be  generalized 
to  that  of  aperture  admittance  by  including  a  radiation  conductance  term. 

If  this  is  not  done,  infinities  in  the  aperture  power  transmitted  may  occur. 

The  result  that  the  transmission  cross  section  of  a  small  aperture 
resonated  by  a  capacitive  load  is  independent  of  the  size  and  shape  of  the 
aperture  is  not  unexpected,  since  similar  results  have  been  obtained  in 
other  transmission  problems.  For  example,  the  transmission  width  of  a 
narrow  infinitely  long  slot  in  a  thick  conductor  at  resonance  is  X/tt, 


J 
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regardless  of  the  actual  slot  width,  as  shown  in  Chapter  VI,  Of  course, 
this  is  for  an  ideal  loss-free  problem.  Conductor  losses  can  significantly 
decrease  the  aperture  transmission  cross  section  for  apertures  in  actual 
conductors . 

The  simplest  way  to  resonate  a  small  aperture  is  to  connect  a 
capacitor  directly  across  the  center  of  the  aperture,  as  was  done  in  the 
example.  This  is  an  analogue  to  the  problem  of  connecting  an  inductor 
in  series  with  a  short  dipole  scatterer.  When  the  aperture  is  in  a  shield, 
we  obviously  do  not  want  to  resonate  the  aperture.  However,  when  a  broad 
spectrum  of  frequencies  is  present,  such  as  in  an  electromagnetic  impulse, 
an  aperture-body  system  will  possibly  resonate  at  some  frequencies  of 
interest.  This  possibility  should  be  taken  into  account  in  the  engineering 
analysis  of  such  a  system. 
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is  provided  to  ESP  Program  Prices  IPOs )  and  other  ESP 
elements.  The  principal  technical  mission  areas  are 
communications ,  electromagnetic  guidance  and  control,  sur- 
veillance  of  ground  and  aerospace  objects,  intelligence  data 
collection  and  handling,  information  system  technology, 
ionospheric  propagation,  solid  state  sciences,  microvoave 
physics  and  electronic  reliability,  maintainability  and 
compatibility. 


